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ABSTRACT 


This design handbook presents data and methods for use in analyzing the dynamical 
performance of a rotor supported in fluid film type bearings. The deta ere given 


in form of design charts covering: 


lL. The dynemic bearirg coefficients 
2. The onset of hydrodynamic instability 
3. The critical speeds and the undalance response of « simple 


rotor-bearing system 


The design charts are in dimensionless form and cover a wide range of bearing 
 aameaie and operating crnditions. 
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’ Several different types of f:uid film bearings are considered together with 


y various forms of lubricatio:: 
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The Plain Cylindrical Gas Bearing, L/D @ 1 «...-.....-.---.-2-..-.-- 
The Plain Cylindrical Gas Bearing, L/D @ 2 ----<-----2c--2-e0n-0-0- 

Critical Speeds and Unbalence Response 

a- Critical Speeds E; 
Plain Cylindrical Bearings, L/D = 0.5, Laminar Film  ..........-..- 
Plain Cylindrical Bearings, L/D = 1, Leminar Film  <-<s<-sesnnennne 
4-Axial Groove Bearings, L/D = 0.5, Laminar Film. pe abuleaTscwees S 
4-Axial Groove Bearings, L/D = 1, Laminar Film  --- -~--.----~-..... 
Elliptical Bearings, L/D = 0.5, == 0.25,. Laminar Film «<---......-- 
Elliptical Bearings, L/D = 0.5, m= 0.5, Laminar Film  «.--....-..- 
Elliptical Bearings, L/D = 1, m= 0.25, Laminar Pilm ------.......- 
Elliptical Bearings, L/D = 1, m= 0.5, Laminar Film  --------..... 


i 

| 
The Plain Cylindrical and the 100 Degree Partial Bearing, Re = 0--- { 
The Plain Cylindrical and the 100 Degree Partial Bearing, Re = 1663 - 
The Plain Cylindrical and th. 100 Degree Partial Bearing, Re = 2378 - 
The Plain Cylindrical and the 100 Degree Partial Bearing, Re = 3326 - | 
The Plain Cylindrical and the 100 Degree Partiel Bearing, Re = $820 - ) 
The Plain Cylindrical and the 100 Degree Partial Bearing, Re = 8315 - 
The Plain Cylindrical and the 100 Degree Partial Bearing, Re = 00 
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D-9 6 Shoe Tilting Pad Bearings, Load between Pads, Leminar Film ......... | 
D-10 5 Shoe Tilting Pad Bearings, Load between Pade, Leminar Film ......... 
D-11 4 Shoe Tilting Pad Bearings, Load between Pads, Laminar Film ......... | 
D-12 6 Shoe Tilting Pad Bearings, Load on Ped, Laminar Pilm ........-ascce | 
D-13 5 Shoe Tilting Pad Bearings, Load on Pad, Laminar Film  ......-2escce 
D-14 4 Shoe Tilting Pad Bearings, '--d on Pad, Laminar Film ...........c08 
D-15 100 Degree Partial Bearings, Alent Film, Re © 1663 -....c..cecncce! 
D-16 100 Degree Partial Bearings, Turbulent Film, Re = 5820 a 
D-17 100 Degree Partial Bearings, Turbulent Film, Re © 13,304 ........... 
D-18 Plain Cylindrical Bearings, Turbulent Film, Re = 1663 testes cath 
D-19 Plain Cylindrical Bearings, Turbulent Film, Re = 5820 Spaciuewecessdul 
D-20 Plain Cylindrical Bearings, Turbulent Film, Re © 13,304 ou vasiea 


| 
b- Amplitude at the Criticat Speeds due to Unbalence 


an a 


i 
D-21 Plain Cylindrical Bearings, L/D = 0.5, Laminar Film ------.--.-. aoe 
D-22 Plain Cylindrical Bearings, L/D = 1, Laminar Film ..-......... ram 
D-23 4-Axial Groove Bearings, L/D = 0.5, Laminar Fila See 
D-24 4-Axial Groove Bearings, L/D = 1, Laminar File wee neceecrercses crows, 
D-25 Elliptical Bearings, L/D = 0.5, m= 0.25, Laminar File ce 
D-26 Elliptical Bearings, L/D = 0.5, m= 0.5, Laminar Film ..............0 
D-27 Elliptical Bearings, L/D = 1, m= 0.25, Laminar Film decentedstedesee 
D-28 Elliptical Bearings, L/D = 1, m= 0.5, Laminar Film encseeeccosenene] 
| 


| 


D-29 6 Shoe Tilting Pad Bearings, Load between Pads, Laminar Film ......./' 
D-30 5 Shoe Tilting Pad Bearings, Load between Pads, Laminar Film  .......! 
D-31 4 Shoe Tilting Pad Bearings, Load between Pads, Laminar Film oo econ) 
D-32 6 Shoe Tilting Pad Bearings, Load on Pad, Laminar Film  .............4! 
D-33 5 Shoe Tilting Pad Bearings, Load on Pad, Laminar Film eee ee 
D-34 4 Shoe Tilting Pad Bearings, Load on Pad, Laminar Film ........... wa! 
D-35 100 Degree Partial Bearings, Turbulent Film, Re © 1663 ...........-..) 
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D-36 100 Degree Partial Bearings, Turbulent Film, Re = 5820 setecse ne cee ee! 


D-37 100 Degree Partial Bearings, Turbulent Film, Re = 13,3046 seccece ewes 
D-38 Plain Cylindrical Bearings, Turbulent Film, Re = 1663 aa 
D-39 Plain Cylindrical Bearings, Turbulent Film, Re = 5820 mesancnecsnnnad 
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D-40 Plain Cylindrical Bearings, Turbulent Film, Re = 13,304 sacs 
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SyMEOLS 
a Orifice radius, inch . 
: 
: B Damping coefficient for hydrostatic bearing, lbe.sec./in 
: Cc Radial bearing clearance, i.e. the difference between the 


radius of curvature of the bearing surface and the journal 
radius, inch 


Clearance of pivot circle in the tilting pad bearing, 
see Page 65, inch 


es 


C. Searing damping coefficients, lbe.sec/in 
D Journal diameter, inch 
d , Feeder hole dianuter, inch 


e Eccentricity of journal center with respect to beering center, i 


| 
| 
| 
} 
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I Mass moment of igertia of pivoted shoe around Longitudinal 
axis, lbe.in.sec 


Tt Transverse mass moment of inertia of rotor around CG, lbe.in.sec 
ub Polar mass moment of inertia of rotor, lbe.in.sec*: 


Bearing stiffness, lbs/in 


Rotor stiffness, lbs/in 


ee BO | 
a On 


Bearing support stiffness, lbs/in 
Bearing spring coefficients, lbs/in 


Bearing length, inch 


moon on 
* 


Distance between admission planes in a hydrostatic bearing, incl 


. ‘ } 

C Ly 2 L-L)> combined length of hydrostatic bearing outside admissiog 

planes, inch : 
C 8 Rotor span between bearings, inch 


Rotor mass per bearing, see Page 202, ibs sec’/ in 


Mo rit = r/e’, Equivalent mass of pivoted shoe, lbs .sec?/in 


m 


Ellipticity of elliptical bearing, see Fig. 8. 
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Rotor speed, RPS 


Total number of feeder holes in hydrostatic bearing (on the charts 


N is used instead of §) 

Ambient pressure, psis 

Supply pressure for hydrostatic bearing, peta 
2 SD, journal radius, inch 

= £ONDC/u, Reynolds number 


Gas constant, in*/sec?.°R 


= wRL(B, 2 Sommerfeld number 

Total temperature, p 

Volume of feeder hole below ortfice, in? 
Bearing load, lbs. 

Rotor emplitude, see Fig. on Page 231, inch 
Journal center coordinates, inch 


Components of journal center velocity, in/sec 


* a + aio for first critical speed and static unbslance } in/ 


= a - Gio for second critical speed and dynamic unbalance 
Influence coefficients for the rotor, see Page 231, in/lbe 
Eccentricity of unbalanced mass, see Figure on Page 23), in. 
= a2/dc, Inherent compensation factor for hydrostatic besring 
= e/C, Eccentricity ratio 
= SERUM By 2 Compressibility number 
a 
6 Na? 
= 3 » Restrictor coefficient 
PC ¥i+d8 
ry 
Lubricant viscosity, lbe.sec/in” 
Whirl frequency, radians/sec 


= 1, for first critical speed and static unbalance 


= (distance between rotor masses) /{,see Figure on Page 231, 
xiii 
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if g - x for hydrostatic bearing 
e e eo Rotor flexibility parameter 


=) 


Bearing attitude angle 


= 2N, Angular speed of rotor, radians/sec 
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7) First or second critical speed, radians/sec 
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¢ The trend towards high-speed turbomachinery haa accentuated the fact that fluid 
film bearings possess stiffness and damping properties which have an appreciable 
effect on the dynamic behavior of the rotor-bearing system. Thus, it has become 
recognized that the static properties (i.e. load carrying capacity, flow, power 
lose, temperature rise, etc) are only some considerations in bearing design. In 
fact, the dynamic characteristics will often govern the particular bearing 
selection. It has been shown that the dynamic properties will affect the threshold 
of instability, critical speeds, unbalance response, noise attenuation, etc. 
Each bearing type (e.g. cylindrical, axial groove, elliptical, pivoted shoe, 
etc), lubrication system (.e.g. self-acting, hybrid), fluid property (e.g. 
compressible, incompressible, and two-phase), and regime of operation (.e.g i 
laminar, transition, and turbulent) will influence the dynamic performance of the 
system to a differing degree. Therefore, in choosing the bearing type for a 


ca 


Us cae 


specific opersting condition and environment, one must consider both the dynamic 


The calculation of the dynamic properties of bearings has been extremely diffi- 
cult and laborious and only the advent of high-speed, high storage computers 

has made it practical. These solutions, however have lagged considerably behind 
steady-state analyses. While at present there are several bocks and numerous | 


| 
; 
and the static aspects. | 
{ 
} 
| 


reports and papers available that treat the steady-state performance and design 
of bearings, only « few papers are available on dynamic properties. This 
handbook partially” fills this gap. The text is presented in such a form that 
as solutions for other bearing types are made available, they can be readily 
added to this handbook. The complexity of the analysis and the numerous parameters . 
also introduce the problem of generalization of the results so as to make them 
broadly applicable. In this text the generalization has been accomplished 
through the use of dimensionless parameters. Decoupling of the bearings from 
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ee 


the rotor system has also been important in order to make the results applicable 


to a variety of rotor configurations. 


—_ 


0 * The scope of this program permitted only a limited coverage of bearing types, | 
' dimensions (L/D, m, arc length, etc) and ranges of operating parameters i 
(] (Re, A, P./P, etc). 
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While the trend to high speed machinery has required the establishment of the 
dynamic properties of bearings, the presence of more severe enviromental re- 
quirements (.e.g extreme temperatures and pressures, rediation, etc), has slso 
focused on this need. High speed aircraft, space power systems (Brayton and 
Rankine dynamic conversion machinery), cryogenic refrigeration machinery, guid- 
ance and control equipment, centrifuges, etc. are just a few examples of new 
machinery requirements. Many of these use process fluid lubrication such as 
liquid metal,:.gas or vapor. The use of these fluids in high speed machinery is 
relrtively recent; therefore, uxperience is limited, and in many cases unfavorable. 
Meuy of the troubles can be attributed to poor design or trial and error procedure. 
At the eame time it should be recognized that the analysis with these process 
fluids is considerably more complex because of the introduction of nonlinearity 
and et least one additional parameter (Reynolds number, compressibility number, 
quality) has to be considered. To ensure satisfactory operation and to prevent 


failures of the machinery it is essential, as part of the actual design procedure, — 


to perform an investigation of the dynamics of the rotor-bearing sys tem. Such 
an investigation must include the calculation of the system critical speeds, a 
calculated prediction of the maximum vibratory amplitude and a check on the. 
stability of the rotor. Through the calculations it can be determined if the 
rotor has adequate stiffness for its weight and if the chosen bearing types and 
the selected bearing dimensions are adequate for the operations] requirements. 
It is the purpose of this handbook to set forth the required design methods and 
to give as much design data as is currently available in easy to use form. 


The data are given in form of design charts which employ dimensionless perameters. 
Thereby, the charts become generally valid and apply to a wide range of operating 
conditions and design configurations. The charts cover three types of data: 


1. Dynamic coefficients for several bearing types and bearing geometries 
(spring and damping coefficients). The data can be used directly 
together with existing computer programe to calculate the unbalance } 
response and the stability of a general rotor-bearing system. 


2. The regions of stable and unstable operation of a rotor bearing system. 


3. The effect of fluid film bearings on the critical speeds and the 
unbalance response of a rotor-bearing system. 
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The first chapter in the handbook contains a brief review of the most important 
aspects of rotor-bearing dynamics from a design point of view. This is followed 
by a chapter in which the usual procedure for arriving et a design of a rotor- 
bearing system is outlined. The various steps in the design procedure are then 
treated separately in the subsequent chapters which at the seme time give a de- 
tailed explanation for using the accompanying design charts. 


That part of the design procedure which deals with the stea:. state performance 
of the system, or the strength and mechanical integrity of the rotor and its 
support structure, has only been included to the extent that it makes the use 

of che presented data more efficient. A more complete coverage of the associated 
design problems can be found in standard referance books and numerous reports. 


Most of the given design data have been computed specifically for the present 
handbook. This includes many cf the charts for the éynamic bearing coefficients, 
almost all of the stability data and all of the charts for the critical speeds 
and the unbalance response. The balance of the data for the dynamic bearing 
coefficients has been compiled from several sources given in the reference list. 


The handbook treats those bearing types which ere classified as fluid film 
bearings. Degign data for ball bearings are given in an accompanying report. 
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‘* The principal part in any piece of rotating machinery is the rotor itself. 


Ite 


function is to generate or transmit power. It consists of a shaft on which may 
be mounted such components as turbine wheels, compressor wheels, gears, the roter 
of an electric motor or generator, etc. The shaft is frequently made of an inte- 
gral piece of material, although the diameter may vary along the length of the 
shaft, and the components are either integral with, shrunk on or mechanically 
fastened to the shaft. The shaft is supported in bearings with separate the 
rotating surface of the shaft from the stationary supporte. The bearings may be 
either fluid film bearings, using ofl, other liquids, gas, or even two phase 
fluids as a lubricant, or they may be rolling element beerings like ball bearings. 


The rotor is never completely rigid and in many applications it is actually 
quite: flexible. Hence, the rotor vibrates if it is subjected to excitational 
forces. These forces may be aerodynamic, magnetic, mechanical, etc. by far 
the most common of such excitations are the unbalance forces present in the rotor. 
They generally derive from the manufacturing processes where it is not possible 
to ensure that the center of gravity of the rotor masses coincide exactly with 
the geometric symmetry axis of the rotor. The unbalance may be distributed 
throughout the rotor and varies in both magnitude and angular position along 

the length of the rotor. As the rotor is brought up in speed the centrifugal 
forces due to the unbalance cause the rotor to deflect and the bent rotor whiris 
around its neutrel axis synchronous with the speed of rotation. This phenomenon 
is known as synchronous whirl. It should be noted that synchronous whirl is 
actually not a vibration of the rotor in the normal sense of the word. The de- 
flected shape uf the rotor itself remains unchanged during the motion and it 


is only when the whirl amplitude is measured in any fixed direction (for instance, 


by a probe) that the motion appears as a vibration. However, from the point of 


view of describing synchronous whirl the analogy to a vibratory motion is very 
strong. This analogy becomes evident when considering the phenomenon of the 


eritical speed of the rotor. For the purpose of illustration assume the rotor to 


consist of a single heavy disc mounted on a flexible shaft whose weight is yery 
small compared to that of the disc. Let the mass of the disc be M in 


and let the stiffness of the rotor be K ibs. (i. e. x. is . 
4 


r in 


i [i the shaft). Furthermore, let the mass center of gravity of the disc be eccentric 

by a distance & inch with respect to the axis of the shaft and let the angular 

‘ c speed of the rotor be w Eadhtne The shaft center at the disc will be whirling 
“ in an orbit with radius x inch. such that the center of gravity of the disc will 

be whirling with a radius of (xta). Hence, the shaft pulle the disc inverds with 

a force of: Ko x lbs. and the centrifugal force pulls outwards with « force of: 


My’ (ta) lbs. For equilibrium the two forces must be equal: 


a 


determined as the force which causes a one inch deflection at the center of 
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Neutral Axis 


Center of Gravity 4 
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It is immediately seen that when the rotor speed is such that: 


ono, IE 


from which the whirl radius of the shaft center at the disc can be computed. 
then the denominator becomes zervu and x becomes infinitely large. The corres- 


ponding speed o, is called the critical speed of the rotor. It may further be 


mm 


noted that for wm the radius x is positive, f.e., x and e are in the same 
direction, such that the center of gravity of the disc is outside the orbit 
of the center of the shaft. For speeds above the critical speed x becomes 


negative and the center of gravity is inside the shaft center orbit. Hence, 


in gaing through the critical speed there is an inversion of the center of 
gravity or,in the terminology of vibrations thephase angle between the ampli- 


a 


tude and the exciting force (i.e. the unbalance) changes from zero to 180 degrees: 
As the rotor speed becomes very large, then x = -a and the center of gravity ; 
of the disccoincides with the center of the shaft orbit and is therefore stationary. © 
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The theoretically infinitely large amplitude at the critical speed does not 
occur in practice since there is always some damping present in the systen. 
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. Cant film forcing separation of the journal and the sleeve. Under steady-state 
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However, any inherent material damping in the shaft itself does not contribute 
to limit the amplitude since the shape of the deflected rotur does not change 
during the whirl motion. Instead, the major source of damping is in most cases 
the journal bearings and without this or a similar source it would be very difficel: 
to pass through the critical speed. It is, therefoze, evident that the bearings — 
play a major role in the dynamics of the rotor. 


The most common journal bearing type is the fluid film bearing in whirh the rotor 
journal and the besring sleeve are separated by a thin film of lubricant. The 
lubricent may be oil, water, liquid metal, or a two phase fluid gas. In the case 


of self-acting bearings the geometry is such that the rotation of the journal cause 
the bearing to act as a viscous pump whereby pressures are generated in the lubri- | 


aS. Helin ae atest ag: 


conditions the total pressure force equals the static load on the bearing. How- 


, ever, if the center of the rotating journal is in motion, as for instance, during 


{ 
synchronous whirl, additional pressures are set up in the lubricant film which 
act as dynamic forces ou the journal in addition to the static force. The dynamic | 
force depends on both the amplitude and the velocity of the journal center motion 
but in cortrast to conventional mechanical forces the dynamic force does not have | 
the same direction as the imposed amplitude. Resolving the dynamic force into 

two components, F. and Fy? along fixed coordirate axes in the bearing, say an 


“y-axis and a veaata:, and likewise resolving the journal center motion into x = 


y displacements, the dynamic force components may be expressed by: 


Fy Kyat ~ Cyy* - yyy > Cyyy 


where x = 3 (x-velocity) and y = x (y-velocity) . To illustrate: 
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Bearing Center 


journal ctr 
position 
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Strictly speaking the given equations are only exact for infinitesimal small 
anplitudes but in practice they prove to be vaiid even for amplitudes as large 


eoccsacel’ 


as a third of the bearing clearance. The 8 coefficients: ne c 
are called spring and damping coefficients such that, for instance, K,, is the 


spring coefficient yielding a restraining force in the x-direction due to a 
displacement in the y-direction. These coefficients are calculated from lubri- 


ae ne 


cation theory and are properties of the particular bearing. They depend on 
the bearing configuration, the lubricant properties and, more important, they bi 
to a given steady-state journal center position. Consequently, they change wi th, 
the speed of the rotor. 


To illustrate the effect of the bearing on the dynamics of the rotor assume 
for simplicity that the bearing may be represented by a single spring co- 
efficient K and a single damping coefficient C instead of the actual 8 co- 

efficients. Returning to the single-dise rotor discussed above the modified 


model now becomes: 


} 
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The critical speed is approximately determined by: 


/ (Kk. : : 
i yet? » 5 cx 


which is less than the critical speed previously calculated with rigid bearings 


i 
where: @, “ff . Hence, the flexibility of the fluid film lowers the critical | 


M 
speed. It is not uncommon for the reduction to amount to 30 to 40 percent and 


cases can be cited where the reduction is even greater, e.g. with soft gas bear- 
ings. Although ignored in the above equation the damping may cause the reducti 

to be less than stated. This is readily seen when considering a bearing with so 
much damping that the bearing acts like it was rigid in which case the critical 

speed becomes equal to st - This effect is frequentiy of importance in oil 
lubricated bearings. On“ the other hand, in gas bearing ‘ applications the rotor — 
is very stiff compared to the gas film in which case the damping causes a slight ! 
reduction in the critical speed. 


When the flexibility of the bearing film is taken into account it is found that 
for whole rotor partakes in the synchronous whirl motions although the sise of 
the whirl orbit changes along the length of the rotor. Furthermore, the presence: 
of damping in the bearing film limits the amplitude: in passing through the - J 
critical speed. The damping is frequently sufficiently large that the critical | 
speed may hardly be noticed and even with instrumentation it may be difficult to 
identify and measure the critical speed with better accuracy than + 10 percent 

in rotor speed. This, however, is also due to the fact that there are 8 bearing | 
coefficients instead of just two. With two coefficients as assumed above the 

stiffness in the vertical and the horizontal direction are the same and likewise 
for the damping coefficients. Hence, due to the symmetry of the system the wht. ef 
orbit is circular. When there are 8 bearing coefficients the system is unsymme- \ 
trical and the whirl orbit becomes elliptical. In essence, this implies that the: 
stiffness and damping coefficients in two mutually perpendicular directicns are | 
dissimilar. Therefore, there will be one critical speed based on the smallest 
stiffness and another critical speed based on the largest stiffness. Theoretical) 
then, there are two resonant amplitude peaks but due to damping the peaks overlap 
sufficiently that in practice they appear as one broad peak. Consequently, it ig 


difficult to define by measurements the exact rotor speed at which the critical 
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speed is encountered. Neither is it possible to define the critical speed by 
means of a change in the phase angle between unbalance and amplitude. The odes 
phase angle value depends on the direction in which the amplitude fs measured. 


Figures 1 to 4 show typical examples of synchronous whirl orbits and rotor ampli 
The amplitudes are measured by capacitance probes which are mounted inside or jue 
adjacent to the bearing. The probes are 90 degrees apart and give the rotor aapl 
tudes as sine-vaves as show by the bottom picture in Fig. 1. By simultaneously } 
displaying the signals from the two probes on an oscilloscope a picture of the 
orbit is obtained on the screen as shown by the top picture in Fig. 1. Hence, 
probes can be thought of as coordinate axis for measuring the rotor motion and 
way be labeled x-probe and y-probe in analogy to x-axie and y-axis. To illustrate 
aseume that the x-signal and the y-signal a-e 90 degrees out of phase so that —_ 
may be written: 


x= a cost y=b cos (wt-90° = b sinct 
i ' (a = peak x-amplitude, b = peak y-amplitude, » = frequency, Tadians, | 
Combining the signals yields: 
xiv 
+ e l 
Pi v2 


i.e. an ellipse with major axis: a(along the x-axis) and minor ams: b(along the : 
y-axis). If the signals are either in phase or 180 degrees out of phase the | 
whirl orbit becomes a straight line. For phuse angle values between 0 and 180 
degrees the orbit is an ellipse whose orientation with respect to the axes depend | 
onthe particular value of the phase angie and on the megnitude of the peak aupli- | 
tudes. This is clearly illustrated by Figs. 2,3 and 4. Figure 2 (and also Fig.l 
is obtained with a rotor supported iz. hybrid gas bearings externally pressurized | 
with rotation)and shows the synchronous whirl orbit, due to unbalance, as a func 
of rotor speed. The amplitude reaches a peak at approximately 18,000 RPM which 1 
either the first or the second crf{tical speed (the two critical speeds are the 
rigid body modes and they are tov close together to seperate). The increase in 
amplitude beyond 25,000 RPM is caused by a pedestal resonance. Apart from the | 
change in size it may also be noted that the orientation of the elliptical orbit 
changes with speed. Figures 3 & 4 present measured synchronous whirl orbits sents 
with a 100 degree, Mlicone ofl lubricated, partial bearing and with a fully 


turbulent lubricant film (Reynolds number = 5620 and 8314, respectively) The, 
measurements are compared to theoretical reaults which are calculated on the bests j 
of theoretically derived spring end damping coefficients representing the bearing ‘ 
film. The overall agreement is very good. Each orbit is shown in the 2-y-coordt- | 
nate system corresponding to the two probes and is obtained for a particular value 
of the Sommerfeld Mumber us @? where: is the lubricant viscosity, Rie the ro 
speed in rps, R is the journal radius, C is the radial clearance and P is the bear- 
ing loading per unit projected bearing area. Hence, the Sommerfeld number gives 
a direct indication of the rotor speed. 


Whereas Fige. 1 to 4 show synchronous whirl with elliptical orbits Fig. 3 showe 
the onset of hydrodynamic instability in a gas lubricated journal beering. The 
orbit has two lobes and is generated by a simple superposition of synchronous 
whirl and half-frequency whirl. The small lobe “corresponds” to the synchronous 
whirl and disappears as the amplitude builds up. If the speed is increased 
further the orbit increases rapidly in size and eventually the bearing will fail 
due to contact between the journal and the bearing. With gas bearings an in- 
stability failure can be very serious and may cause extensive damage to both the . 
rotor and the bearings and even with liquid or oil lubricated bearings it ie very 
difficult to operate a rotor beyond the onset of inetability. The hydrodynamic 
instability is known as either fractional frequency whirl ("half-frequency whirl") | 
or resonant ship ("oil whip") The instability is inherent in the very principle { 
by which lubrication works, To illustrate, consider for simplicity a vertical 
rotor supported in plain journal bearing with length L, radius R, radial clearance | 


C and an angular rotor speed of w Sather. 
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Let the journal ride eccentric tn the. bearing clearance with on eccentricity e. 
sider a flow balance for one half of the lubricant film (shaded in the figure). 
and assume that no pressure is present in the fila: 


flow in= 41 nw (or) ia 


3 
flow out @ 4 LRo (C-e) dn 


Since the flow into the film volume is greater than the flow going out the file 
will normally develop a pressure such that the inflow is reduced and the outflow | 
ie increased with side leakage taking place. It is by this mechanism that a | 
hydrodynamic bearing is able to generate its load carrying capacity. However, | 
assume that the journal center is not stationary but whirls with an instantaneous 
angular speed of v Eadiens around the bearing center. The journal center velocity 
is then: 3°” An. and the considered film volume increases with time at the rate of: | 
2LRev ex’ In this case the flow balance can be written: 


(flow in) = (flow out) + (rate of chads of volume) 


or 
Aino (cre) = +L Ro (C-e) + 2Rev : 
which is satisfied if: | | 
Ps , 
ves 2° \ 


When the above flow balance is satisfied no pressure is developed in the bearing fi 
Hence, when the whirl frequency equals half of the rotational speed the bearing 
looses its load carrying capacity and it becomes unstable. 


| See 


The purpose of the above illustration is to show that whirl may cause instability 
and that the instability whirl frequency has an inherent tendency to equal a 

half of the rotor speed. However, as already mentioned the illustration is only 
valid for a vertical rotor or a lightly loaded bearing where the rotor whirls 
around the bearing center and the whirl orbit is a circle. In general the bear- , 
ing is under static load. Consequently the center of the whirl orbit is dis- \ 


ll 


placed from the bearing center end the orbit becomes elliptical. This hes the | 
effect of reducing the instability whirl frequency sithough the reduction ie only. 
noticeable when the static load and eccentricity is high (for a sufficiently high 
eccentricity, the instability may actually be eliminated in the oil lubricated bear- | 
ing). Furthermore, the simple concept of instability as being caused by a complete 
loss of load carrying capacity is no longer valid for a loaded bearing. Rather, 
the onset of instability can be viewed as being caused by the bearing loosing its 
damping ability. A brief explanation of this phenomenon will be given in the 
following. 


It can be shown that, for the purpose of describing the onset cf instability, . 

@ journal bearing can be represented by a single effective spring coefficient and 
a eingle, effective damping coefficient. These coefficients can be calculated 
from the 8 dynamic bearing cocfficients as shown im Ref. 16. However, the two 
coefficients derived in this manner depend on the whirl frequency such that the 
effective damping coefficient is negative for small frequencies and becomes 
positive for higher frequencies. The frequency at which the damping becomes sero 
shall be called the instability frequency. The ratio between the instability 
frequency and the rotor speed equals approximately 0.5 in most cases except at 
high eccentricity ratios (e.g. at low speeds) where the ratio is usually smaller. 
Therefore, the instability frequency can be found as @ function of speed as shown 
by the corresponding curve in the figure below. Yor frequencies leas than the 
instability frequency (i.e. in the region below the curve) the effective demping 
is negative, and for frequencies greater than the instability frequency the damp- 
ing is positive. The effective spring coefficient of the bearing together with 
the flexibility of the rotor determines the resonant frequencies of the rotor- 
bearing system. Since the bearing stiffness is a functicn of speed the resonant 
frequencies become speed dependent. The lowest of these resonant frequencies ia 
shown by the curve labeled: "System Resonant Fruquency" in the figure below. 
This curve intersects the curve for the instability frequency at a speed denoted 
as: "Instability Threshold Speed": 
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System Resonant 


Rotor Speed 
st Critical Instability e 
Speed Threshold 
Speed 


Assume that the rotor is subjected to a small disturbance. It will then tend to 
vibrate at its lowest resonant frequency. However, if the rotor is running below 
the instability threshold speed the bearings provide positive damping and the 
vibration dies out. As the speed is increased the damping available for the 
vibration diminishes until it becomes equal to zero at the instability threshold 
speed. Attempting to increase the speed beyond the threshold speed causes the 
damping to become negative such that any imposed disturbance is amplified and 

the system io unstable. 


From the figure it is seen that if the rotor mass is increased or the shaft is 
made more flexible the system's natural frequency is lowered, whereby the inter- 
section between the two frequency curves shifts to the left and the threshold 
speed is reduced. Conversely, if external damping is present in the system (for 
instance, in the supports) the curve of the film instability frequency is lowered 
thereby raising the threshold speed. The latter characteristic suggests means 

by which an otherwise unstable system may be stabilized and also explains why 
some machinery, notably rotary pumps, may operate stable well above the theore- 
tical threshold speed. The damping capacity of the liquid passing through the 
impeller acts as external damping to the system and stabilises the rotor. 
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Turning agein to the phenomenon of critical speadsthe preceding discussion has — 
only mentioned the lowest critical speed (the f‘rst critical speed). Since the ; 
rotor is an elastic body there are, of course, infinitely many critical speeds but 
only the lowest two or three critical speeds are of importance in practice. Each 
critical speed has associated with it a mode shape which is simply the shape of 

the deflected rotor at the critical speed. To illustrate the interaction between 
the rotor and its bearings at the critical speeds consider a simple model con- 
sisting of a uniform shaft supported in bearings with equal stiffness, denoted 


' eat 


If the bearings are rigid = «) the shaft behaves like a simply supported 
beam with mode shapes: 


lst Critical Spoed 2nd Critical Speed 3rd Critical Speed 


Rigid Bearings, Flexible Rotor ; 


Next, assume the bearings to be very flexible in comparison with the rotor. Then 
the critical speeds will be lowered and the rotor will behave like a rigid body 
at the first and the second critical speed with mode shapes*) : 


ie = ae 


lst Critical Speed 2nd Critical Speed 
(Translatory Whirl) (Conical Whirl) 


Flexible Bearings, Rigid Rotor | 
(Rigid Body Critical Speeds) 


tile, satin cect wrens ME as 


* Strictly speaking, it is only true that the lst critical speed is the translatory | 
whirl mode and the second critical speed is the conical whirl mode when the trans- 
verse radius of gyration is less than 70 percent of the rotor span between bear- | 
ings. In practice this condition is usually satisfied. 
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It should be remembered that the rotor does not vibrate but whirls. At the 

firet rigid body critical speed the whirling rotor describes a cylindrical 
surface in space and the mode is called: translatory whirl. At the second rigid 
body critical speed the described surface is a cone with its agsx between the 
bearings and the mode ia called: conical whirl. Both of the two rigid body 
critical speeds are proportional to /X. Hence, if the critical speeds were 
plotted as a function of bearing stiffness on a log-log graph the rigid body 
critical speeds will appear as straight lines with « slope of 0.5 and the flexible 
rotor modes will appear as horisontal lines: 


Bearing Stiffness K (log scale 


There are only two rigid body modes. Even if the bearings are very flexible 
the third critical speed, therefore, involves bending of the rotor itself and 


it is the free-free mode: 


3rd Critical Speed 
(free-free mode) 


Flexible Bearings 
For a rotor supported in gas bearings the bearing stiffness is generally suffi- 
ciently small that the two lowest critical speeds are the rigid body modes. 
Actually, at this time it is recommended tn require this condition since gas 
bearings do not possess enough camping to effectively control the resonant 
15 
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amplitude of a flexible rotor. Hence, in current design practice a gas bearing sh, 
supported rotor is quite rigid and the third critical speed is normelly outside 
the operating speed range. ; oe) 


* 


7 . 
The ahowe discussion summarizes the basic elements of rotor dynamics from the 

point of view of design application. The three major subjects are: a) unbalance 
response and synchronous whirl, b) critical speeds and, c) hydrodynamic instability.. 
To this may be added rotor balancing which, howewer, will be discussed in another 
report (Ref.1). Other considerations are in most cases of secondary importance 
although they may of course be very significant in certain special pieces of 
machinery (shock loading, non-synchronous excitation, random vibration, noise 
transmission, etc). The most frequently encountered of the additional considera- 
tions sre treated in Ref. 1 to the extent informetion is available. 
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OUTLINE OF THE DESIGN PROCEDURE FOR A ROTOR-BEARIMG SYSTEM ' 


A successful rotor design for high performance rotating machinery applications 
requires the integration of many technical considerations, one of which is rotor- 
dearing dynamics. Compromises must be made to balance the often conflicting 
requirements of materials properties, stress, heat transfer, aerodynamics, and 
electrical design in addition to rotor-bearing dynamics. Careful attention must 
also be given to the supporting structures in which the bearings are mounted. 
These structures must maintain the aligament of the bearings with the shaft aad 
with each other and maintain clearances around rotating members. At the same 
time they must provide proper structural stiffness for the support of the rotor 
and its dynamic loads. The design should be such that these functions are ac- 
complished with minimum change during load and temperature transients. i 
The remainder of this section is devoted to the rotor-bearing dynamics aspect , 
of rotor design. It should be remembered, however, that this is only one of many 
considerations that must come into play during the evolvement of the rotor design. 
At an early stage of the design it is necessary to take into account the critical 
speeds and stability of the rotor-bearing system, and the methods for balancing 
tue rotor and minimizing rotor amplitude. The following chart illustrates the 
relationship of rotor-bearing dynamics to other design considerations in the: 


-- 
oe 


design procedure. ie a 


Rotor specification 


Specified rotor function 


Configuration of rotor and supporting structure 


ho eee — — Dimensions, stiffness and weight of rotor 
| f | es Selection of journal and thrust bearing type 


a 
| Sub jects 
| treated 
, Calculation of unbalance response in this 
> —-— <@#-—jBalancing specifications and dynamic bearing loads manual. 
L Stress, Heat transfer, aerodynamic, electrical and materials considerations 
[Final design| 
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The “boxes” give the sequeatial steps of the design procedure while the broken 
lines indicate “feed-back" paths depending the outcome of the various calculations. 


The particular machinery application dictates the function of the rotor. Typical 
examples are: turbo-generator, turbo-compressor (gas turbines in general), 
turbo-compressor alternator, turbo-alternator pump, wotor-compressor or motor- 
pump. There are also multi-rotor systems where two or more rotors are coupled 
together as in conventional turbine-generator sets, connected with a flexible or 
a splined coupling, or gear shafting where the gears themselves provide the con- 


' nection. between the individual shafts. The required output of the machinery 


determines the optimum number of turbine or compressor stages and wheels and 
their ’size,. or the size of the rotor and stator for an, electrical motor or genera- 
tor. On this. basis a preliminary lay-out of the rotor can be made incorporating 
such other considerations as the strength of the rotor, thermal insulation, heat 
shielding, space for seals, space for inlet and exhaust ducts, and assembly 
sequence. The approximate weight of the rotor and the corresponding bearing 
reactions can be computed from the preliminary layout. In addition, the rotor 
span, the journal diameters and the space available for the bearings can be 
determined. 


Early in the design process, attention must be given to the stiffness and the 
mass distribution of the rotor. Stiffness is of particular concern if the rotor 
is built-up of several shrink-fitted sleeves, collars and disks compressed | 
together by means of either a nut on the shaft or a long tie bolt. Ideally, the 
built-up rotor is as stiff as an integral rotor if all the parts have exact 
dimensions and the compression between parts induced by the pre-loading is sut-~ 
ficiently high that it wili never be exceeded by any tensile stresses caused by | 
bending_of the rotor. However, this is never achieved in practice and a reduction 
in stiffness is always encountered. It can only be minimized ty having as few 
separate componénts as possible, by requiring tight tolerances and by specifying 
a high pre-load. An’ additional disadvantage of the built-up rotor is the possi- 
bility of eacountering instability induced by dry. friction and sliding betwee: 
the rotor components. Such an instability sets in just above the first critical 


speed. os, 
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The mass distribution of the rotor is important because of its strong influence 
on critical speed. For example, the rotor frequently extends beyond the bear | 
te carry a tuchine wheel, a compressor wheel or a similar relatively heavy mass.; 
This is normally referred to as a rotor with overhung mass. Such an arrangveent | 
causes several problems. The primary reason is that the shaft sections near and 
through the bearings are usually the most flexibie part of the rotor. Thies is | 
especially the case in high speed appiications where the size of the journal dia 
meter is kept to a minimum in order to reduce the friction power loss in the 
bearings. Hence, the overhung masses are supported on rilatively flexible | 
“cantilevers" thereby causing a significant lowering of the critical. speeds. The 
flexible shaft sections at the ends aleo make it more difficult to balance the 
rotor. Normally the only accessible balancing planes are at the end of the rotot 
and if heavy correction weights are required, the rotor will bend and make it | 
virtually impossible to obtain ao satisfactory balance. 
Another effect of overhung mass is to shift the shaft nodes towards the bearings, 
The resulting decrease in whirl amplitude at the bearings reduces the capacity | 
of the beuring to provide damping. Consequently, the whirl amplitude of the rotq 
between the bearings will get very large. This problem is especially severe at 
the third critical speed. 
nee | 
The shape of the overhung mass is also important in the way it influences the 
so-called gyroscopic effect. The dynamic forces on relatively thin disks like 


a turbine wheel tend to exert a stiffening effect on the rotor at high speed 


and raise the critical speed. Relatively long cylindrical masses like an 

alternator rotor have an opposite effect and lower the critical speed. As will | 
be shown in subsequent chapters, the factor which determines which type of | 
behavior will occur is the relative magnitude of the polar and transverse 


moments of inertia >f the mass. 


To arrive at a final rotor design requires detailed calculations and the incorpor 
tion of the bearing characteristics. The choice of bearing type and dimensions 
igs largely governed by the properties of the available lubricant, the cpereting 
environments, the speed of the rotor, the load imposed on the bearings and the 
dynamic performance of the bearings, including stability. A more detailed 
discussion is given in the next chapter which also describes the various bearing 


types for which data are available in this manual. 
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Once the preliminary lay-out of the rotor and its beerings is established, the 
spring and damping coefficients for the bearings can be determined from the 
design charts, (see Pigs. B-1 to 3-128), and « calculation of the critical speeds 
can be performed, usually by oceans of a caaputer program. If the operating speed 
ie too close to any of the critical speeds, the rotor or the bearing design suet 
be modified. To assist in this evaluation, the design charts, Fige. D-l to D-20 
can be used to determine the effect of bearing type and bearing dimensions on the 


critical speeds. 


Knowing the dynamic bearing coefficients also makes it possible to check the. ; 
stability of the rotor-bearing system. For this purpose a computer program may ; 
be employed or the design charts, Figs. C-1 to C-17, can be vsed. This is dis- 


cussed in detail in a later chapter. 


Finally, the magnitude of the rotor whirl amplitude and the corresponding 
dynamic bearing loads due to a rotor unbalance should be estimated. Thereby the 
sensitivity of the rotor to the location of the unbalance can be evaluated and 

the required balance level can be established. The calculations are normally 
performed by a computer program but an overall estimate can be obtained from 
Figs. D-21 to D-40 in which the effect of both static and dynamic unbalance is 


presented. 


The subsequent chapters are devoted to a detailed discussion of each of the y 
steps in the outlined design procedure and to a description of the presented 


design charts. 
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The wide variety of machinery epplications has produced numerous bearing types 
to suit the different requirements. A representative selection of such bearing 
types are treated in the present report as summarized in Fig. 6. A brief des- 
cription of eact type is given below: 


‘a. The Plain Cylindrical Bearing (See Fig. 7) 1s the simplest of all journal | 


bearing types. Its performance characteristics are well established and ex- 
tensive design information is available. Hence, the cylindrical bearing plays | 
an important role in the study and analysis of rotor-bearing systems. Its | 
practical application, on che other hand, is generally limited to the gas lubri-' 
cated bearing. To make the plain cylindrical bearing practical for oil lubri- 
cation or other liquid lubricants it is necessary to provide the bearing with 
grooves or holes through which to supply the lubricant. Sometimes a single, 
circumferential groove in the middle of the bearing is used, and in other cases : 
one or several axial grooves are provided. { 
b. The 4-Axial Groove Bearing (See Fig. 7) is one of the commonly used oil lubri 
cated bearing types. It is a cylindrical bearing provided with four oil supply 
grooves. The oil is supplied at a nominal gage pressure which ensures an adequa, 
oil flow and thereby some cooling of the bearing. Occasionally the grooved bear‘ 
ing has been used in gas bearing applications in order to improve the bearings 
stability but this case is not considered here. 


c 11 r (See Fig. 8) is another typical example of an oil 
lubricated bearing. It is especially used in gear and turbine applications. As: 
a bearing type it is classified as a lobed bearing in contrast to a grooved bear. 
Whereas the grooved bearing consists of a number of partial arcs with a common 
center the lobed bearing is made up of partial arcs whose centers do not coin- 
cide. Thus, the elliptical bearing consists of two partial arcs where the 
bottom arc has its center a distance m: C above the bearing center and vice 
versa for the top arc (see Fig. 8). Here C is the bearing clearance and m is 
called the ellipticity. The value of m is between 0 and l. By this arrangement 
the bearing is preloaded which means that the journal center eccentricity 
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with respect to the loaded arc is increased and never becomes zero. Thereby 
the bearing is stiffened and its stability is sonewhet improved. 
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a. The Partial Bearing (See Fig. 8) is in general not a practical bearing ia 
itself. However, having available partial bearing data allows combining the 
partial arcs into a variety of grooved and lobed bearing configurations. The 
present report treats 4 partial bearing configurations: 50, 60, 80 and 100 degree 
arc, (the centrally loaded partial bearing). The lubricant is in all cases a 
liquid and the bearing film is laminar except for the 100 degree arc where tur- 
bulent flow is also considered. 


a, The Tilting Pad Bearing (See Fig.9) is widely used in high speed applications 


i 
{ 


where hydrodynamic instability and misalignment are problems. The tilting pad bear- 


ing consists of a number of shoes, each shoe actually being partial bearing, 
which are mounted on pivots. Hence, the shoes are free to adjust and to follow 
the motions of the journal. In this way the tilting pad bearing ensures in- 
herent stability as long as the mass inertias of the shoes do not interfere with 


the adjustment ability of the bearing. In the present report the 3,4, 5 and 6-shoe 


tilting pad bearing are considered. The load direction may either pase between 
the two bottom shoes ("load between pads") or it may pass through the pivot of 
the bottom shoe ("load on pad"), (see Fig. 9). The lubricant is incompressible 


(liquid) and the lubricant film is laminar. For the 3 and 4-shoe bearing, turbulent | 


flow is also considered. 


- £. The Hydrostatic Bearing (see Fig. 10) finds application where the bearing load 


is very high, when a very stiff bearing is required or where the rotor speed is 
too low to allow the bearing to operate hydrodynamically. The lubricant is 
pressure fed to the bearing through flow restricted feeding holes and, thus, 

the bearing requires equipment to furnish the flow such as a pump, & compressor 
or an accumulator. In the present case the gas lubricated hydrostatic bearing is 


considered. 


&:. The Ball Bearing offers the advantage of compactness, simple lubrication re- 
quirements and relatively high speed capability. However, its operational life 
is limited. Since the ball bearing as a bearing type is entirely different from 
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the fluid film bearing types the design dats and performance data for the bell 
bearing are given in a separate report. 


Selection of Bearing Type 

At a very early stage of the design it is necessary to choose the bearing type. 
Unfortunately, the choice is governed by many factors and i¢ is dangerous to 
give any general rules on which to base the selection. Most frequently the 
bearing type is chosen on the basis of steady state performance; Icad carrying 
capacity, lubricant flow requirements, friction power loss, etc. but in many 
applications it is also necessary to consider the dynamical per formance; dynamic 
load, stiffness, damping and stability. Since the present report is concerned 
with rotor-bearing dynamics it is not possible within this scope to also treat 


the rather large subject of bearing performance and design for steady atate 
operation. Reference shall instead be given to the published literature among 
which References 3 & 4 may serve as representative examples. However, a brief 
summary will be given of the major considerations affecting the selection 

of the bearing type: ' 4 

1. The Lubricant | 
Examples of lubricants are: oils, gases, liquid metals,water, vapour, etc. The 
choice of lubricant is dictated by what can be made available, the required 
auxiliary equipment, and the operating conditions as discussed later. The lubri- 
cant must be chemical stable in the operating temperature range, and if the : 
environment is radioactive conventional lubricants are affected. From a purely 
lubrication point of view the most important lubricant characteristic is its 
viscosity and its boundary lubrication ability. For oils and most liquids the 
viscosity is strongly dependent on temperature and it decreases with increasing 
temperature. For gases the viscosity increases with temperature but the effect 
is not nearly as pronounced. The density-to-viscosity ratio is sometimes of 
importance especially with low kinematic viscosity fluids such as liquid metals, 
and the higher the ratio the more is the bearing likely to operate in the turbu- 


lent regime. 


2. Rotor Speed 


From a functional point of view there are three most coamonly used bearing 
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types: the hydrodynamic (self-acting) bearing, the hybrid (externally pressarived . 
with rotation) bearing, and rotiiig element trarings > (The first tee types =re 
fluid film bearings where the rotor surface is seperated from the bearing surface 
by a thin film of lubricant. The hydrodynamic bearing works ss a viscous pap 
and its load carrying capacity is, therefore, proportional to the lubricant 
viscosity and the rotor speed. The hydrostatic bearing, on the other hand, 
carries its load by means of the pressure at which the lubricant is supplied. 
Thus, for very low speeds or where the load is very high the hydrostatic bearing 
should be used (alternatively a rolling element bearing) . Wot moderately high 
speeds the liquid lubricated hydrodynamic bearing canbe, used, the latter bearing 
having the advantage at higher speeds due to its lower friction. For very high 
speeds the gas lubricated bearing should be applied. 


2._Bearing Load 


The bearing's load carrying capacity should be considered simul taneously with 

the rotor speed as discussed above. In fluid film bearings the maximum load 
carrying capacity is determined by specifying the minimum allowable filmthickness 
taking into consideration the type of machinery application, possible dynamic loads 
and the size of possible dirt particles. Wit gas bearings and lubricants with 
low boundary lubrication ability there are certain limits imposed on the bearing 
load because of the dry friction and wear taking place during starting and 
stopping before the self-acting film has formed. If the load is ‘too light, on 
the other hand, the bearing may be susceptible to hydrodynamic instability. 


4, Friction Power Lose 
The magnitude of the power loss in the bearing is frequently e major factor in 
deciding which bearing type to use. For fluid film bearings operating in the 
laminar regime the power loss is approximately given by: 


~ 3, 2 
friction power loss = ear eT HP 
6600-C- Vl-« 


where: lubricant viscosity, lbs.sec/ in? 


rotor speed, RPS 
journal diameter, inch 
bearing length, inch 


radial clearance, inch 


a arwg a F 
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eccentricity ratio = (journal center eccentricity) /(radial clearance) 
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For eccentricity ratios ebove .5 to .6 the actual power loss becomes a little 
higher. It fa sean that the power loss increases with the square of the speed 
and the cube of the diameter. Thus, for high speed applications the use of an 
oil as the lubricant would neceseitate a very small journal diameter. Thereby 
the rotor may become too flexible and instead it is necessary to choose a low 
viscosity lubricant like gas or, in some cases, liquid metal. In the latter 

case the bearing film frequently becomes turbulent with a corresponding increase | 


in the power loss. = | 


a. Operating Temperature 


Ig the bearing operates at an elevated temperature a conventional hydrocarbon 
oil bearing is no longer feasible and instead a special oil or a gas bearing 

(in special cases a liquid metal beering) aust be used. For very low temperatures 
oils again are not available and the gas bearing may be applied. Apart from — 
affecting the choice of lubricant the operating temperature also affects the 
selection of the bearing and shaft materials and configuration. Thus, elevated 
temperatures are in general accompanied by appreciable thermal gradients with 
corresponding distortions of the bearing geometry. It is, therefore, important 
to select a bearing type which is not too sensitive to thermal distortions or 
which can be provided with means to ensure a uniform temperature distribution. 


§. Dynamic Load , 5? Sy 
In addition to the static load a bearing is also subjected to some dynamic load, 
for instance, unbalance forces or shock forces. Actually, the magnitude of the 
dynamic forces cannot be established prior to the design of the rotor-bearing 
system since they depend on the rotor mass, the rotor flexibility and the dynamic 
properties of the bearing. Even so, they should be estimated and taken into 
account in selecting the bearing for the application. 


i, Bearing Stiffness and Damping 


The bearing stiffness greatly influences the critical speeds of the rotor-bearing 
system. The critical speeds, in turn, are of major importance in determining 

the operating speed range of the rotor. As a general rule: the higher the number 
of critical speeds in the speed range, the more difficult becomes the operation 
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of the rotor and the more difficult it is to balance the rotor. Furthermore, 
if the unit is not a variable speed machine the design speed should be removed 
from the closest critical speed by at least 10 to 20 per cent. Knowing the 
critical speeds, however, is only sufficient in establishing those speeds at 
which a relatively high rotor whirl amplitude can be expected. To evaluate the 
amplitude it is also necessary to consider the bearing damping. To assist in 
comparing various bearing types on the basis of their influence on the critical 
speeds and the corresponding whirl amplitudes, several design charts have been 


prepared as discussed later. 


8. Rotor-Bearing Stability 

The necessity of avoiding hydrodynamic instability is frequently a determining 
factor in selecting the bearing type. Thus, in most high speed gas bearing 
applications the tilting pad bearing will be preferred over other bearing types 
because of its very high threshold speed. Design charts for determining the 
onset of instability for several bearing types are included in the present re- 
port and will be discussed later. 


In any given application the consideratiors listed above enter into selecting 
the bearing type. The final choice must be based on a numerical evaluation of 
each of the factors. The present report gives data and methods for evaluating 
the factors connected with the dynamical performance and the stability of the 
rotor-bearing system. Data for the steady-state performance is available else- 
where, for instance Refs. 3 & 4. 


Selection of Bearing Dimensions 


In most cases it is difficult to separate the selection of the bearing type from 


the selection of the bearing dimensions. In dimensioning and designing the bear- 


ing the same factors must be considered as previously discussed for selecting 
the bearing type. Among the factors the dynamical performance will be treated 
in subsequent sections and in the following a brief discussion will be given of 
design data for calculating the bearings load carrying capacity. The discussion 
is concerned with fluid film bearings; ball bearings are treated in a separate 


report. 
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The design data for load carrying capacity ere given in Figs.A-l-to A-14. They 
_are in dimensionless form such that they cover a wide variation in design con- 
ditions. The governing dimensionless parameters are summarized in Fig. 11 

and will be further defined below. 


— CA & 


For liquid lubricants (incompressible lubricants; there are 4 dimensionless 


parameters: 


cm 


Length-to-diameter ratio: L/D | 


[| Sommerfeld Number: $= ute (b? Laminar 

Turbulent | 

tf Eccentricity ratio: e= e/C | 

Reynolds Number: Re = sont 

ir | 

[ where: 
i L - Bearing length, inch 


D - Journal d‘ameter, inch 


4 D, Journal radius, inch 


C7 
) 
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- Radial clearance, inch 


- Journal eccentricity, inch 

- Rotor speed, RPS 

Bearing load, lbs. 

- Lubricant viscosity, Lba.sec/in” ; 


- Lubricant mass density, lbs.sec*/in 


The radial clearance C is always the difference between the radius of curvature 


ae | 


of the bearing surface and the journal radius R. The journal eccentricity e 


is the distance between the journal center and the bearing center. 


a 


For the elliptical bearing an additional parameter enters, namely the ellipticity 
m. It is defined as the distance between the lobe center and the bearing center 
divided by the radial clearance C, see Fig. 8. Similarly, in a tilting pad bear- 


care Ui S| 


ing the center of curvature for the shoes may not coincide with the bearing center 


ne a 


in which case a similar parameter is required as for the elliptical bearing. For 


the tilting pad bearing this parameter is known as the preload factor: (1-c'/Cc). 


This factor is zero in all the design data for load carrying capacity. 
27 
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Figs. A-1 to A-1l give the eccentricity ratio or the miniaam filethickeess as 4 
function of the Sommerfeld number for various bearing types end geometries, all | 

with liquid lubricants. Under specified operating conditions and with known 

bearing dimensions the Scamerfeld number can be computed and the cerresponding 
eccentricity ratio (or minimum filmthickness) can be found from the charts. In 

the cases where the eccentricity ratio is given, the minimm filathickness can [ 


be found from: ee 

minimum film thickness = C - (l-«) [: 
Yor each of the three tilting pad bearings there are 2 charts. The first chart [ 
(Figs. A-6, A-8 and A-10) gives the relationship between eccentricity ratio and ; 
Sommerfeld number. The eccentricity ratio is denoted ag eR to emphasize that [ 
it applies to the complete bearing, not to the individual pads. The second 2 
chart (Figs. A-7, A-9 and A-11) gives the filmthickness over the pivot for the [, 


bottom pads as a function of eR: 


The charts for partial bearings (Figs. A-3 to A-5) may occasionally be useful 

in estimating the eccentricity ratio and minimum filmthickneas for milti-groove [., 
or multi-lobe bearings. For multi-groove bearings a close estimate is obtained : 
by assuming that the loads is carried solely by the bottom segment of the grooved 
bearing. To illustrate, it is seen from Fig. 7 that the bottom segment of the 
present 4-axial groove bearing configuration is 90 degrees. Hence, the eccentricity | 
ratio for the 4-axial groove bearing, already covered byFig. A-l, can also be ob- 

tained with good accuracy from Fig. A-5 which is valid for an 80 degree partial | 
bearing. A similar procedure can be employed for a multi-lobe bearing but only 

at relatively high eccentricity ratios. 


| 
It should be emphasized, however, that the outlined method must be used with great | 
caution and judgement since all the partial bearing data in the present handbc ok 
assume a centrally loaded pad. Furthermore, extensive data for eccentricity ratio 
and minimum filmthickness are already available elsewhere for several multi-groove 
and multi-lobe bearing configurations (see Ref. 3). ‘Thus, the chief purpose of the [ 
presented partial bearing data is not so much to define the steady-state performance 
of the bearing but, rather, to make it possible to determine the eccentricity ratio | 
at which the dynamic bearing coefficients are evaluated. This will be discussed 

| 


later. 
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For the self-acting gas bearing the dimensionless design parameters are: 


Length-to-diameter ratio: L/D 


I 
1 
F 


; Compressibility number: Aes Maul &? 

é a 

i Eccentricity ratio: «= e/C 

: 

; Bearing load parameter: W/P LD . 


where: 
Pr. - Ambient Pressure, psia 


ge 
- 
s 


and the other symbols have been defined above. 


oe of 


Figs. A-12 to A-14 allow determining the eccentricity ratio for the plain 
cylindrical gas bearing when the values of the compressibility number A and 


~~ 


the load parameter are known. : a 


reas | 


For the hydrostatic gas bearing the operating eccentricity ratio can be cal- 


culated with good approximation from: 


Ww m 
e= Xx (e < .5) 


where K is the stiffness in lbs/in as determined from Pigs. B-110 to 8-122. 
The calculation of K will be discussed later. The eccentricity ratio computed 


in this way should nut exceed .4 to .5 in order to avoid lock-up. 


SS 


“\ 
of the minimum filmthickness. The choice is based on experience and is governed 


of cooling of the bearing and the maximum tolerable bearing surface temperature. 
As a rule, high speed ratios (e.g. centrifuges, high speed gears, superchargers, 
etc) will operate with a low eccentricity ratio (€ < .4) while low speed gears, 
turbines, etc. will operate with a relatively high eccentricity- ratio («€ = .4 to 


.8). For further comments, consult Refs. 3 and 4. 


On this basis the given design charts together with the formulae for the 
friction power loss can be used at least to arrive at a preliminary set of 


bearing dimensions. 
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In general it is difficult to give specific rules for the’ smallest allowable value 


by such considerations as the size of possible dirt particles, the desired amount 
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Bearing Ne. 2, Probes: 
Supply Pressure = 200 psig 


Shaft Orbit 


Scale: Vertical Axis - 24C. in./major div. 
Horizontal Axis - 206, in./aajor div. 


Top: horizontal Amplitude (x-probe) 
Bottom: Vertical Amplitude (y-probe) 


Bearing-to-Shaft, 50 su/cm 
Anbient Pressure = 14.7 psia 


Speed = 3C 300 RPM 


Fig. 1 


Synchronous Whirl due to Mechanical Unbalance 
Hydrostatic Gas Bearing 


Whirl Orbit and Corresponding Amplitude Components 
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Fig. 2 Synchronous Whirl due to Mechanical Unbalance 
Hydrostatic Gas Bearing 
Whirl Orbit as a Function of Rotor Speed 
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TYPES OF BEARINGS CONSIDERED 
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@ PLAIN CYLINDRICAL BEARING: LIQUID LUBRICATED- LAMINAR ANO TURBULENT FILM 
GAS LUBRICATED 


Are er aere Re 


b 4-AXIAL GROOVE BEARING LIQUID LUBRICATED-LAMINAR FILM 


c. ELLIPTICAL BEARING LIQUID LUBRICATED -LAMINAR FILM 
ob : 
ae d. PARTIAL BEARING LIQUID LUBRICATED -LAMINAR AND TURBULENT FUM 
ce : | 
aes [. e. TILTING PAD BEARING LIQUID LUBRICATED -LAMINAR AND TURBULENT FILM 
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f. HYDROSTATIC BEARING GAS LUBRICATED 
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7 [ ¢. BALL BEARING 
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| [ Fig. 6 A Schematic Summary of the Bearing Types 


for which Data are given in the Report 
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4-AXIAL GROOVE BEARING 


Pig. 7 Geometry of the Plain Cylindrical Bearing 
. and the 4-Axial Groove Bearing 
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Fig. 9 Geometry of the Tilting Pad Bearing 
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EQUIVALENT DYNAMICAL MODEL 
OF ROTOR-BEARING- PEDESTAL SYSTEM 


Fig. 12 Dynamical Model of Rotor~Bearing Systen 


REPRESENTATION OF THE BEARING FILM 
UNDER DYNAMIC CONDITIONS 
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Fig. 13. Dynamical Model of the Bearing Fluid Fila 
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Figure A-7 The 6 Shoe Tilting Pad Bearing, Load between Pads . 
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Figure A-9 The 5 Shoe Tilting Pad Bearing, Load between Pads 
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Figure A-12 The Plain Cylindrical Gas Bearing 
Load vs. Compressibility Number and Eccentricity Ratiu 
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Figure A-13 The Plain Cylindrical Gas Bearing 
Load va. Compressibility Number and Eccentricity Ratio 
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Figure A-14 The Plain Cylindrical Gas Bearing 
Load vs. Compressibility Numter and Eccentricity Ratio 
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SPRING AND DAMPING COEFFICIENTS FOR FLUID FJ lef JOURNAL BEARINGS 


The journal bearings supporting a rotor may be either of two general types: 
rolling element bearings or fluid film bear.ngs. In both cases the bearings 
possess flexibility and damping and these properties influence the dynamic 
characteristics of the rotor. Thus, in order to perform a design analysis of 

the rotur-bearing system it is necessary to know the values of the bearing stiff- 
ness and damping. For this purpose the present report gives desigr: data tor e 
large number of fluid film bearing ‘types. The data are in tom of graphs, as 
Fig. B-l to B-124 which will be discussed in detail below. Rolling element 
bearings (ball bearings) are treated in a separate report (Ref. 2). 


General Concept of Spring and Damping Coefficients 


In all fluid film bearings the journal is separated from the bearin surface 

by a thin film of lubricant. The lubricant may be an oil, eT Fave or 
a ges; and the film may operate in the laminar or the turbulent flow regime. 
In hydrodynamic (self-acting) bearings the load carrying capacity is due to pressures 
generated in the fluid film and this in turn requires the journal to operate 
eccentric within the bearing or, in other words, the larger the eccentricity 
the higher the load carrying capacity, everything else kept the same (including 
the speed). This may be illustrated by considering one of the simplest bear- 
ing types namely the plain cylindrical bearing. It consists of @ cylindrical 
sleeve with radius R and with length L. The difference between the. bearing 

” radius and the journal radius is the radial clearance C. Thus, the distance 

between the bearing center and the journal center, denoted as the eccentricity e, 
can vary between zero and C. For convenience the eccentricity is normalized 
with respect to the clearance, the ratio being known as the eccentricity ratio 

| ¢ = e/C. Hence, € varies between 0 and 1. Assume the rotor to be running with 
a speed of N rps and apply a load W to the bearing, starting from zero and 


increasing steadily. ‘then the journal center will describe a Locus: 
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8. DYNAMIC BEARING \ 
COEFFICIENTS 


Rearing Center 


Journal Center 


Journal Center Locus 0 


x 
The position of the journal center with respect to the bearing center is given 


by the eccentricity ratio ¢ and the attitude angle @ (see figure above) such 


that @ goes from approximately 90° to 0° as ¢ goes from 0 to l. The vertical 
displacement in the same direction as the load becomes ecos¢. 


Introduce an x-y-coordinate system with the x-axis vertically downwards whereby 
x = ecos® = Cecos¢ and y = esin@ = Cesiné. Hence, the static stiffness is given 


by: 
ae 
Static stiffness dx ” Cd(ecos¢) 


f.e., it is the slope of the load curve plotted as a function of the vertical 
displacement. Since the load curve is not linear with displacement the stiff- 
ness changes with eccentricity ratio too: 
Static 
Stiffness 
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It is evident that the static stiffness alone is not sufficient to describe the 

total etiffnese characteristics of the bearing. From the sbove journal center : 
locus it is seen that even under static load the journal displacement is not in 

the same direction as the appiied load. The displacement has both a vertical [. 
component x and a horizontal component y. This is in contrast to conventional , 
mechanical systems where the elastic derormation usually is in the same direction | 

as the force. Hence, it is necessary to split up the static stiffness in two 

parts, one part related to the x-displacement and one part related to the y- [ 
displacement. Assume the rotor to be operating at a given speed under a static ; 
load W and with journal center coordinates x and y. Increase the load by an [ 
infinitesimal small amount A whereby the journal center coordinates are changed 

by Ox and Ay, respectively. However, to separate the effects of Ax and Ay, the | 

change should be performed in two steps: first increase x by Ax such that W is ‘ 
increased by MW, and next, increase y by Oy whereby W is futher increased by [, 


Ou. In total: 


= My + a, 


or oa, ou, 
Ad = at ax + a « by : 


For very small changes the ratios become derivatives. Since the derivatives are 
taken with respect to one coordinate at a time, they are partial derivatives: [| 


fa” Oy 

—_l+ —2v ie 

ae iy * by Lj 
' Therefore: [| 


g 


o He ax + HE ay 


This may be shown to be a separation into components of the previously discussed 


em — = 
ean 


atatic stiffness as follows: 


aa vy Oa fa. dy 
dx Ax zt oi y dx 
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For convenience the partial derivatives are designated by symbols: 
wey 
xx Cx xy y 


They are called spring coefficients. There are, however, more than two coeffi- 


cients as will te discussed. 


Under static conditicas thc static bearing load W is opposed by an equal and 
opposite force W,. representing the integrated pressures in the fluid film. Wy 
is taken to be positive in the same direction as x such that: 
WH) static % -W 
If, however, the journal center is outside the earlier shown journal center locus 
the fluid film pressures generate an additional horizontal force component W 
which under static conditions forces the journal back to its equilibrium position 
where W. = 0. This means that the bearing has stiffness not only in the vertical 
direction (x-direction) but also in the horizontal direction (y-direction). 
Hence, in calculating the partial derivatives above where x is first increasd by 
4x, not only is W increasd by a) (more accurately it is W,. which is increased 
by o,) but Wy is changed by ou. Similarly is Wy changed by fe, in increasing 
y by Oy. Therefore, in total: 
Lat Ce ow ow | 
(ND seatic™ “ * = ans 5 : fg => / dx * 5# iii 


si sa ow ow 
CHD geatte RB IE RT et oy 


The minus signs are due to the convention of Wand Wy being positive in the same 
direction as x and y, respectively. It is convenient to write these equations 


as; 
(ai) static = Ky ea Ry oy 


(ty) static = Kyx oe Ky” ey 
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where: x «= oe ees OW 
, xx xy —— t 
oy nf : 
os = K ao 
yx yy dy | 
These four coefficients are the spring coefficients of the bearing. The first 


index gives the direction of the force, the second index gives the direction of 
the displacement. It is seen that a displacement in, say the x-direction, gives | 
rise to a force not only in the x-direction, but also in the y-direction. 


Under dynamic conditions the journal center no longer occupies a fixed position. 

It has a velocity with components x and y (x # as and y «= & vhere t ie time). 
Such a velocity induces additional pressures in the fluid film resulting in an 

increase of the fluid film force components We and ¥,: This effect is called [i 
the squeeze film effect because the forces are generated by the lubricants re- 

sistance to being "squeezed" out of the bearing film. In complete analogy to | 
the four spring coefficients the squeeze film forces can be represented by four 

damping coefficients: Ce cy? Ca and C yy’ In total, the changes in the [ 
fluid film force components due to a motion of the journal center are given by: 


on se ae ax - Ce x= Ky Ly - Cy -y 


au, ae Key > hx - Coe x Ky - Ay - Cyy °¥y | 


The coefficients where the indices are the same (i.e. Bo Ky and C . are | 
called the direct-coupling terms (force and displacement are in the same direction) 
whereas the coefficients where the indices are different (i.e. Ky xy Kye and ci. , 
are called the cross-coupling terms. Since the cross-coupling terms do not possess the 
necessary symmetry it is not possible by a simt.c rotation of the x-y-coordinate | 


system to eliminate the cross-coupling terms. 


coupling spring coefficient is not a pure stiffness, it also contains damping 
properties. This mey be illustrated by considering the energy dissipated during 


Furthermore, the cross-coupling terms are not “pure” coefficients, i.e. a crose- 
the journal center motion. Assume the journal center to whirl in an orbit ( 


with angular frequency w such that: 


i 


€ } € . + ‘ 


é e ’ ' e + t ! : } 


} 4 é i } | } . ; ! : : j 3 } oy 


Ax © a cos(wt + a) i.e. x @ - Wa sin(wt + *) 


Sy = b sin(wt + o) y= wb cort(wt + *) 


The energy dissipated during one orbit of revolution is: 


y 


tL : 

= \ : «|. er 

Energy dissipated ti ae &x + Ce + K yy + Pb as c y 
u 


Ke K ‘yt#C ylysde 
x y? vo) 


2 2 
- “ - - . @ -0 ~— 
= Ka we x b aoe a.b. sin( * Awcxy + ur) ab - cos( re y "xy y I 
Hence, the cross-coupling spring coefficients may contribute to the energy dissi- 


pation which is usually considered a prope<ty of a damping coefficient. 


The 8 dynamic coefficients are calculated from lubrication theory. They are computed 
as the previously indicated partial derivatives of the camponents uf the fluid film 
force. They depend not only og the bearing type and the bearing geometry but also on 
the steady state journal center position. ‘Thus, to each operating point on the journal 
center locus (see the figure, page 59) corresponds different values of the 8 coeffi- 
cients. Therefore, when the eccentricity ratio is known the dynamic coefficienta are 
also defined. The eccentricity ratio in turn is a function of several dimensionless 


parameters: 


for liquid lubricated, hydrodynamic bearings: 


Length-to-diameter ratio : L/D 


Sommerfeld Number: S= BNDL (&? laminar 
W Cc urbulent 
Reynolds Number : Re = sONDC/p 


(for laminar film: Re = 0) 
; t 


for gas lubricated, hydrodynamic bearings: 


Length to-diameter ratio: L/D 
Compressibility Number : A» Saul &? 

a 
Bearing load parameter: w/P LD 
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for gas lubricated, hydrostatic bearings: 
Length tu-diameter ratio: L/D 
A ai 6uN VR T | 


Restrictor coefficient: 


Supply Pressure Ratio: Pi/P. 


The symbols are defined on Page 27 and in the nomenclature list. 


Use of the Design Charts 


The design data for spring and damping coefficients are given in Figs. B-l to 
B-124 for a wide selection of bearing types. The bearing types considered are 


listed in Fig.6. 


a, Liquid Lubricant, Laminar Film 


The design data for the bearing types with an incompressible lubricant and a 
laminar film are given in Figs. B-l to B-55. The charts give the dynamic co- 
efficients as a function of the Sommerfeld number S. The coefficients are in 


dimensionless form: 


CK CK CK CK 
—th —aAY. —YA ‘a. | 
W W W W 
CuwC CC Cac Cw ] 
—— KK — YL —— YK —LY | 
W W W W 
where: | 


C - Radial clearance, inch 

W - Bearing loed, lbs. 

® - 2nN. Angular speed, radians/sec | 
N - Rotor speed, RPS . | 
K_,K_,K__,K__ - Spring coefficients, lbs/in 

C_,C _,C _,C_ = Damping coefficients, Lbs.sec/in | 


When the bearing dimensions, the load and the viscosity of the lubricant are known | 
the Sommerfeld number can be computed for several values of the speed over the 


operating range. Entering the charts with these calculated values fhe correspond- | 
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ing values of the coefficients can be computed. Thereby the coefficients 
become determined over the speed range and can be used in an undalance 


response calculation, a stability analysis or a critical speed calcuation. 


For the plain cylindrical, the 4-axial groove and the partial bearing the 
only other parameter besides the Sommerfeld number is the L/D = ratio. For 
the elliptical bearing an additional parameter is the ellipticity m which 
gives the distance between the lobe center and the bearing center normalized 
with respect to the radial clearance. Thus, m= 0 corresponds to a plain 
cylindrical bearing whereas m = 1 would mean solid contact between the 
journal and the bearing. In the present design charts (Figs. B-9 to B-16) 
two values of the ellipticity are considered: m= 0.25 and m= 0.5. 


In analogy with the elliptical bearing the tilting pad bearing may not have 


the center cf curvature of the pads coincide with the center of the bearing. 
If the journal radius is R and the radial clearance is C then the radius 

of curvature of the pads is (R+C). Let the circle with center in the 
bearing center and touching the surface of the pads at the pivots have 

' + the radius (R+C').« If C' = C the pad centers coincide with the bearing 
center and if C' = 0 the journal is always in solid contact with the pads. 
The parameter (1 - C'/C) is called the preload factor and it may vary from 
0 to 1 like the ellipticity m. 


In Figs. B-37 to B-52 the preload factor is zero. Furthermore, Figs. B-37 to 
B-48 assume that the mass inertia of the shoes is very small and can be 
ignored. Under this assumption the 4 cross-coupling terms of the dynamic 
coefficients vanish and oniy K_, ,C__ and C__ remain. In addition, for 
xx" yy” Xx yy 
the 4 shoe bearing with the load direction passing between the two bottom 
pads K = Ky and Cs = Cyy because of symmetry. The effect of pad mass 
inertia and the corresponding values of the cross-coupling terms are pre- 
sented in Figs. B-49 to B-52 for the 4-shoe tilting pad bearing. The mass 


inertia is expressed in terms of a dimensionless parameter: 


—— 
for?) 


; 3 i 3 . ‘ 7 X : ‘ ’ ’ ‘ ' : i : : . , a . ae . : i ‘ . . é e ae 


di 


where: 


M = r/R? 
I = mass moment of inertia of a pad around a lgngitudinal 
axis passing through the pivot, lbs.in.sec 


The other symbols are defined on Page 27. 


In addition to the dynamic coefficients the design charts for the tilting pad 
bearing also show a curve labeled “critical mass". For each Soumerfeld number 
this curve gives the value of the above dimensionless masa inertia parameter 
at which the phase angle between the journal motion and the pad motion becomes 
90 degrees. If the actual pad inertia is more than approximately half the 
critical value the cross-coupling terms become significant. The critical mass 
parameter also gives some measure of the bearings stability characteristics 


as discussed later. 


The three last charts, Figs. B-53 to B-55, apply to a vertical rotor and show | 
the effect of the preload factor. Since the load W is equal to zero for a 
vertical rotor it is necessary to redefine the dimensionless form of the dynamic 
coefficients. The revised form is defined on the charts. 


The values of the dynamic coefficients for other multi-greove bearings than the. 
4-axial groove bearing and for other multi-lobe bearings than the elliptical | 
bearing can be estimated approximately from the partial bearing coefficients 
given in Figs. B-17 to B-36. This is best illustrated by firet giving the 
exact and general method for calculating the dynamic coefficients of a bearing 
made up of several segments, say n segments (i.e.for the 4-axial groove bearing! 
n= 4, and for the elliptical bearing: n= 2). Consider the i'th seguent, whic 
is actually a partial bearing in itself, and let this segment be located an ang! 


vy, with respect to the vertical load line for the complete bearing. vy is meaat 
“in the direction of the shaft rotation from the top of the bearing to the radia! 


line bisecting the segment. Next, it is necessary to know ths steady-state loci 
of the journal center within the complete bearing for the given load direction. 


. This information can be obtained from Refs. 3 or 4. Once the journas center l¢ 


is available the eccentricity ratio and attitude angle for each segment can be 
established from simple geometrical relationships which depend on the particula 


‘ 
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bearing configuration under investigation. Knowing the eccentricity ratio and 
the attitude angle for the i'th segment allows computing the 8 dynamic coeffi- 

: eibtarers ri 
cients for the segment Roa ere (es from the analysis deacri 
Combining the n segments results in the 8 dynamic coefficients fq 


the complete bearing as follows: 


previously. 


n 
Kx "5 [x .cos¥+ i sin’¥ - ( K cos¥ sid J t 


n 
Rey 12 [x,,co8"v- er sine + (K.7 Ky) cont siny J ; 


n 
ee "5 [x,,co#'v- Ky sin’y + ( Ky cosy sid ] t 


K = : [r cos V+ K_ sin’¥ + (K_+K cos¥ sind | 
yy ih yy xx xy sie. 1 
The equations for we? Oe ys arn and wy are the same by simply replacing 
K with wC in the four equetions above. 


_It is evident that the segments in general are not centrally loaded. Hence, i 
the partial bearing data given in this handbook are tc be used for calculating 
‘the dynamic coefficients of a multi-groove or multi-lobe cearing certain assum 
tions are required. First, for the particular eteady-state opevating conditia 
the eccentricity ratio and the attitude angle for each segment must be known 
(e.g. frov kefs. 3 or 4). Thereby the load line for each segment is also know 
and the angle vy should be measured to this line instead of the bisecting line 
Next, the load line dividea the segment into two arcs which 
Multiply the smallest arc by two and use this arc lengt 
as the effective partial bearing arc for the segment. If the load line is out 
With the obtained values of eccentric 


as used before. 


normally are unequal. 


side the segment disregard the segment. 
racio and effective arc length enter the appropriate chart from Fig. A-3 to A-' 
read off the Sommerfeld number and use the corresponding chart from Figs. Bel? 
to B-34 to determine the 8 dynamic coefficients for the segment. Repeat this 

procedure for each segment and, finally, compute the overall dynamic coefficie 


for the complete bearing from the equations given above. 
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The condition of the flow in the bearing film, whether laminar or turbulent, 
is determined by the Reynolds number: 


Re = “AIRC 
rv 


where: 
~ journal diameter, inch 


- radial clearance, inch 

Totor speed, RPS 

~ lubricant viscosity, lbs.sec/ia” 

~ mass density of lubricant, lbe.sec?/in’ 


oF BOA 
a 


When Re < 1l000thefilm is laminar. For higher Reynolds number the film becomes 


turbulent. 


The dynamic bearing coefficients are given for three different bearing types 
operating in the turbulent flow regime: the plain cylindrical bearing, the 

100 degree partial bearing, the 4-shoe tilting pad bearing and the 3-shoe tilt- 
ing pad bearing. The data are summarized in design charts, see Figs. 8-56 to 
B-93. Apart from the Reynolds number the dimensionless parameter are the same. 
as for the laminar film data. 


Cc. b d 


When the lubricant is compressible the Sommerfeld number is no longer adequate 
as the governing design parameter. Instead, the compressibility number Ais 


used: ‘ 
_ A2uuN (BR 2 
eT sank: 


where: . 
Pr. - Ambient pressure, psia 


and the other symbols are defined on Page 27. In addition, one more parameter 
ie necessary. It can either be the eccentricity ratio « or the dimensionless 
bearing load wW/P LD (note: Sommerfeld number S = A/(12x - w/P LD) ). In the 
present deaign charts for the dynamic bearing coefficients the. chosen parameters 
are Aand e¢€. The relationship between A, ¢ and w/P LD is shown in Figs. A-12 
to A-14, 
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The dynamic coefficients are given in dimensionless forn: | 


— AL 
P LD P_LD P LD P_ LD 
a a a a 


cK cK cK cK 
—m YK _y | 
PLD PLD PLD PLD 
a a a a 
Cue Cw cnc | 
—m th: ee 
i 


Hence, when the bearing dimensions, the load, the rotor speed and the gas visca 
are known the compressibility number, A, and the dimensionless load can he comp 
Entering Figa. A-12 to A-14 the corresponding value of the eccentricity ratio ¢ 
be found. With this value and the value of A the dynamic coefficients can be _ 
determined from Figs. B-94 to B-109. Repeating the calculations for several spi 
allows determining the coefficients over the complete speed range. They can tly 
be used in calculating the unbalance response and the critical speeds of the ro 


bearing systen. 


When the lubricant is a gas the dynamic coefficients are actualiy a function of 
the whirl frequency in addition to depending on Aand ¢. In the present case u 
whirl frequency is set equal to the rotational speed of the journal. This 

that the design charts are only valid for calculations involving synchronous wht 


i.e. unbalance response anu critical speeds. 


a. Gas Lubricant, Hydrostatic Bearing 
Whereas all the previously considered bearing types are self-acting and rely 
on the viscous pump effect (hydrodynamic effect) for generating their load capa 
the hydrostatic gas bearing is pressurized from an external source (a compressed 
or a bottle). Therefore, the presentation of design data for the hydrostatic 
bearing requires a new set of design perameters. These parametcrs are: 

& for single plane admission 
Lengtheto-diameter: & = . , 


a 
D for double plane admission 


Restrictor coefficient: A, ae 3 


Supply pressure ratio: Po/P 
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Yeeder hole volume ratio: NV. /aDLC 


where: 


D - Journal diameter, inch 

R ¢- 4 D, Journal radius, inch 

L - Bearing length, inch 

C - Radial clearance, inch 

Ly - Distance between admission planes, inch 

L, - L-L,, Combined length outside adaission planes, inch 

P. - Ambient pressure, psis 

Ee - Supply pressure, psia 

a - Orifice radius, inch 

d= Feeder hole diameter, inch . as ptt. 

N - Total number of feeder holes (in the design charts N is used 
in the place of R) 

Vv, 7 Volume of feeder hole below orifice, in? 

R - Gas constant, in?/eec*. oR 

T - Total temperature, oR 

uw 060) SC Gas viscosity, Ibeceec/iae 

& - a*/dc, Inherent compensation factor 


(for details, see Fig. 10) 


Instead of 8 coefficients as for a hydrodynamic bearing the hydrostatic bearing 
has only two coefficients: a spring coefficient K and a damping coefficient B. 


They are given in dimensionless form: 


Dimensionless stiffness: a8 . _-_._ 
14% 3¢ (P.- P,) DL 


Dimensionless damping: ui? 


where: 
K © Spring coefficient, lbs/in 


B - Damping coefficient, lbs.sec/in 


The hydrostatic gas journal bearing is a plain cylindrical bearing provided with 
feeder holee for admission of the pressurized gas. The feeder holes may be 


arranged in one or two planes, denoted as single plane or double plane admission. 


70 


eee cE mmc apennereer ronan a ene RETA I Rep RARE ee ote .g 
+ Be PE 
; ; ree 


ee oe | 


| 
The feeder holes are flow restricted by means of an orifice and some restriction 
may also take place in the “curtain area" between the rim of the feeder hole and 
the journal surface. The latter area is equal to xdC whereas the srea of the 
orifice is Kar: The ratio between the two areas is: 5 = a2/dc which, therefore, 
gives an indication of where the major part of the flow restriction takes place. 
The restrictor coefficient A, gives the ratio between the flow resistance of the 
bearing film and the flow resistance of the feeder holes. If the film resistance, 
is very small (A, = 0) the pressure in the gas film becomes equal to ambient and 
the bearing has no stiffness or load carrying capacity. On the other hand, if thy 
film resistance becomes very high or the feeder hole resistance is very small 
(A ~»=) the pressure downstream of-all the feeder holes will cqual the supply 
pressure and again the bearing has no stiffness. There is, therefore, an optimau 


value of the restrictor ratio, approximately given by: 


Ae * .5 to .7 


Under this condition the stiffness is a maximum. 

The restrictor coefficient together with the supply pressure ratio are the con- 
trolling desigr. parameters for a hydrostatic bearing. However, in considering 
the bearing damping a third one enters namely, the feeder hole volume ratio. Th 
parameter gives the ratio between the volume contained in all the feeder holes 
and the volume of the gas film. It affects the time lag between the flow and 
pressure when the journal performs dynamic motions. The value of the parameter ; 
should in general be less than 05 to .1 to minimize the bearings susceptibility 


to pneumatic hammer instability. 


For given bearing dimensions and lubricant p erties the controlling design par 
meters may be computed and Figs. B-110 to B- can be entered to determine the 
stiffness and damping of the hydrostatic bearing. The values given for the dang 
ing coefficient assume that the vibratory frequency is not too large. 


In an actual application the journal rotation superimposes hydrodynamic pressur¢ 
on the hydrostatic pressure. This effect, usually denoted as the hybrid bearing 
effect, has been ignored in the present design charts. The effect is normally | 


small. n 
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Figure B-1 The Plain Cylindrical Bearing, Laminar File 
Dimensionless Dynamic Coefficients 
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Figure B-2 The Plain cylindrical Bearing, Laminar Film 
Dimensionless Dynamic Coefficients 
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Figure B-22 The 50 Degree Pactial Bearing, Centrally Loaded, Laminar Fil 
Dimensionless Dynamic. Coefficients 
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Figure B-23 The 60 Degree Partial Bearing, Centrally Loaded, Laminar Film 
Dimensionless Dynamic Coefficients 


ee eee 


96 


ow 


C 
C 
E 
E 
f 
i 


10 S=SS=S5 Pe ee 33 pete +R RUSK. 3: SS SE EW RERLS 31th: WC BBB Ziti 
e =S=SS=5 BISS3Sz SSSI UE SS Re MB SBSS555 1582 SS 258981152: SSB ISILI 5s MBE is 
EES Bess HH PS reiting 3ISI5355 E323 mi Beue: 
‘ BEE EIBGBRE 
7 Cn oe Gn En ore 4 3 eo ee 
SSS Se SSS SENN ES it 
eee Rota rE 
SSSceriiie uuu BORN BH f 
2 nm 
SSRSSSEISSH REBEL Sa 
(Gh GP GED OE GD SO SHHSIGHEDS Seeeseqsne added 
er oe Th @SBe 00tsn oo cha 
Smmaaacssiiinaiiers nore 
ttt ty BRRSTR. ena 
TT CHT HHH 14 To 
MMMBBR SRST An 
ans bated bed ee ti 


fae SPL 887i wa SBI 2s? 
SRS Tere: Ba TH ss 8s 881: 


Pe Lh rth 
ch sina 


i itt ia 
rn jiinae iW Hot HH i A 
CCC REET nnn Hi i” Ht Hit at 
CCC eC i 
gC COCCI ALERT aa ! vu H wi I 
0 10°! | 10 104 


Figure B-24 The 60 Degree Partial Bearing, Centrally Loaded, Laminar Fiin | 
Dimensionless Dynamic Coefficients 
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Figure B-25 The 60 Degree Partial Bearing, Centrally Loaded, Laminar Film 
Dimensionless Dynamic Coefficients 
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Figure B-26 The 60 Degree Partial Bearing, Centrally Loaded, Laminar Fila 
Dimensionless Dynamic Coefficients 
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Figure B-27 The 60 Degree Partial Bearing, Centrally Loeded, Laminar Fila 
Dimensionless Dynamic Coefficients 
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Figure B-28 The 60 Degree Partial Bearing, Centrally Loaded, Laminar Fils 
Dimensionless Dynamic Coefficients 
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Figure B-29 The 80 Degree Partial Bearing, Centrally Loaded, Laminar FI! 
Dimensionless Dynami.c Coefficients 
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Figure B-30 The 80 Degree Partial Bearing, Centrally Loaded, Laminar 1 
Dimensionless Dynamic Coefficients 
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Figure B-32 The 80 Degree Partial Bearing, Centrally Loaded, Laminar Fil 
Dimensionless Dynamic Coefficients 
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Figure B-33 The 80 Degree Partial Bearing, Centrally Loaded, Laminar Fila 
Dimensionless Dynamic Coefficients 


106 


i 


veo 


‘ 


Co an 


prc 


i 


f 1 Ne pa ee 
t a: : ve a 


é 


¢ . f ter, — 


; ‘ : é ' 


s 


L 


me ee ee ee ee ae er ee ees ism estos: 
PSSST Trt. +. tt ttt er titt tt de Tt Priel — ee 
PP tereriiis st ts tee Serre Seco mee aot me es Saletan a me we 
208 800) eam com a aes st: to Rose cesses: ewe Dawes ern: ae 
06e 02 een He a OSES Ese foclaedy-dndab-beedt Lvndvanel ADT RESHN i> Soe 
Co Se pds mene, 
= 3 saeers Shae As Semen ae oon: 5 Sekt Ladelebed wage: 
ae 298 84: Hit mooaasn 50 name Se one Saneacsssiaes Senses: 
ee eee eee eee 096001 ana am msoarenes simeeases | 
Lidbe ae 


eee tea mer aeet naan 


2a a wows oe 


seecis 


etic 


b+ 
sstecuaeusace Lbdhed : 


2 | ae 
meen 
RS". cht 
Fesiisul ire 
Bnet itt 
4 sae 


Hee 
H 


i 
Ht atl HH 
Ny 


a 
igs 
iz 
siz: : = 32322 
a HE HH Sa 


ay | | scclimauenis: 
smeeesriua J.tLHUN aimee saensuirime seein 
z= #=5 sezaie: 

Oe Ore ae me ae we 
tj HT auoucee 

+t 4 me ees eans 

Ld 0) am oe aan eeee 
am r] dae @uensar 

is vam @easias 

ize mamas’ we t+] Na aasates 
ae a eben 

f imi 

il 


Hh HTH it 
aa sthiga 
Pest 

CES 


be 


Sossscen 


Figure B-34 The 80 Degree Partial Bearing, Centrally Loaded, Laminar I 
Dimensionless Dynamic Coefficients 
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Figure 3-35 The 100 Degree Partial Bearing, Centrally Loaded, Laminar Fila 
Dimensionless Dynamic Coefficients 
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Figure B-36 The 100 Degree Partial Bearing, Centrally Loaded, Laminar Fila 
Dimensionless Dynamic Coefficients 
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Figure B-37 The 6 Shoe Tilting Pad Bearing, L/D = .25, Laminar Film 
Load between Pads 
c'/C = 1, No Pad Inertia 
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Figure B-38 The 6 Shoe Tilting Pad Bearing, L/D = .5, Laminar Film 
Load between Pads 
c'/C # 1, No Pad Inertia 
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Figure B-39 The 6 Shoe Tilting Pad Bearing, L/D @ :.75, Laminar File 
Load between Pads 
C'/C = 1, No Pad Inertia 
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Pigure B-40 The 5 Shoe Tilting Pad Bearing, L/D = .25, Leminar Fila 


Load between Pads 
c'/C = 1, No Pad Inertia 
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Figure B-41 The 5 Shoe Tilting Pad Bearing, L/D # .5, Laminar Film 


Load between Pacis 


c'/C = 1, No Pad Inertia 
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Figure B-42 The 5 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film 
Load between Pads 
c'/C = 1, No Pad Inert‘a 
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Figure B-43 The 4 Shoe Tilting Pad Bearing, L/D = .5, Laminar Film 


Load between Pads 


c'/c = 1, No Pad Inertia 


116 


RO. a leanne 


\ CRIT. MASS | 


: | ae are Fed a 


& il CQ) 
di Pe 
i Cwlxy YH “1 Azle 
ei wl PHO 
Pa Cel eg 
| Na | 


| 
\ 
| 


Figure B-44 The 4 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film 
Load between Pads 
c'/C = 1, No Paj Inertia 
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Figure B-45 The 4 Shoe Tilting Pad Bearing, L/D = 1, Laminar Film 
Load between Fads 
c'/c #1, 


No Pad Inertia 
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Figure B-47 The 5 Shoe Tilting Pad Bearing, L/D = .5, Laminar Film 
Load on Pad 
c'/C = 1, No Pad Inertia 
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Figure B-48 The 4 Shoe Tilting Pad Bearing, L/D = 


Load on Pad 
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c'/C = 1, No Pad Inertia 
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Figure B-49 The 4 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film 
Load between Pads 
c'/C = 1, Effect of Pad Inertia 
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Figure B-50 The 4 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film 
Load between Pads 

c'/C = 1, Effect of Pad Inertia 


| 123 


| 
1. | 
|) Se 710 | an | | 


yet Fao 2] a Bi 
a 
4 | | Be 


oom 
TE th 
BEil 


10" a 
CTT SIT 


“) af 2 46 wl 2 468] 7 
2 
g: HO 


Figure B-5! The 4 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film 
Load between Pads 
c'/C = 1, Effect of Pad Inertia 
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Figure B-53 The 12 Shoe Tilting Pad Bearing, L/D = .25, Laminar Film 
Vertical Rotor, Effect of Preload 
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Figure B-54 The 6 Shoe Tilting Pad Bearing, L/D @ .5, Laminar Film 
Vertical Rotor, Erfect of Preload 


127 


Figure B-55 The 4 Shoe Tilting Pad Bearing, L/D = .75, Laminar Film | 
Vertical Rotor, Effect of Preload : 
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Figure 3-56 The Plain Cylindrical Bearing, Turbulent Film 
Dimensionless Dynamic Coefficients 
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Figure B-57 The Plain Cylindrical Bearing, Turbulent Fila 
Dimensionless Dynamic Coefficients 
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Figure B-77 The 100 Degree Partial Bearing, Centrally Loaded, Turbulent Film 
Dimensionless Dynamic Coefficients 
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The 100 Degree Partial Bearing, 


Figure B-78 
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Figure B-79 _The 100 Degree Partial Bearing, Centrally Loaded, Turbulent Film 
Dimensionless Dynamic Coefficients 
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Figure B-80 The 4 Shoe Tilting Pad Bearing, Turbulent Pils [| 
Load between Pads 
c'/C = 1, No Pad Inertia fi 
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Figure B-82 The 4 Shoe Tilting Pad Bearing, Turbulent Fila 
Load between Pads 
c'/c = 1, No Pad Inertia 
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Figure B-83 The 4 Shoe Tilting Pad Bearing, Turbulent Film 
Load between Pads im 
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Figure B-85 The 4 Shoe Tilting Pad Bearing, Turbulent Fila 
Load between Pads 
c'/C = 1, Ho Pad Inertia 
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Figure B-86 Tha 4 Shoe Tilting Pad Bearing, Turbulent FA lu 
Load between Pads 
c'/C = 1, Mo Pad Inertia 
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Figure B-87 The 3 Shoe Tilting Pad Bearing, Turbulent Fila 
Load between Pads 
c'/C = l, No Pad Inertia 
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: Figure B-88 The 3 Shoe Tilting Pad Bearing, Turbulent Filz 
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Figure B-89 The 3 Shoe Tilting Pad Bearing, Turbulent Film 


Load between Pads 
c'/C = 1, No Pad Inertia 
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Figure B-90 The 3 Shoe Tilting Pad Bearing, Turbulent Fila [, 
Load between Pads 
c'/C = 1, No Pad Inertia 
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Figure B-91 The 3 Shoe Tilting Pad Bearing, 


Load between Pads 
C'/C = 1, No Pad Inertia 


165 


dons hadeth-Lobah 
Ht tote teri am Bea tats) 


ity 


Baars) 


manages Sosa stele 
Heim eR Ett 


HHH 


} Hiei, 


ant ian stir, 
GUTH EOLA: 


«WOR ineart 


O0500: Min Gm Od. WH E80 188s: 


RESETS Me Titty 
tte Seeae8 


Turbulent File 


HNL MARA NIE 


aura: 


rs 
vie eee 
TSP OSER & 
ei ti aa] 
TsIGae es 
Hi 


Mts 380588 Sew eee 
OE ECC ROU St. KERR 
te A Pt ‘2288 Sag 
ACCT 
Es 
% 


1eunes 


pape 
iY 


Ss 
- 
- 
- 
= 
= 

x 


- 
es 


SOSUULE MB SERIA 


sie 
abeaiiie 


! 
; 
* 
é 
‘ 
+ 
‘ 
t 


agettit 
O° adaent 
20050001 


a 


4 
2 
. 
a 
* 
€ 
' 
' 


” 


OS RAUSING AT BAS SHH s HTD. 


t 
' 
1 
f 
} 


ai uu 


PAISittene 
NORE TET 


peas er 


¢rsscome 

revsaaeeam 
ieteee CORR eg wh enteose 
MICCERRO Ma litestegm 
ACCOR RE OS Meesseee 
HC UCCROR Om istteege 
10 0ORE Rw Merteses 


ng, Turbulent Film 


Li] 
88881 Mw Seas aees 


+ 


Hts 
166 


mbaeit 


OReTALH 
c'/c =.1, No Pad Inertia 


OR wages 


SBuit: 
=m 
i 
ww egeery 


a 


PE Ty 


euane 
Figure B-92 The 3 Shoe Tilting Pad Beari 
Load between Pads 


GileeeaP 
aT 
BEER 
Silitiaeaas 


mteenenet um 
abeeer een ae 


ieestesse 
eeeaaean 


7 La 
agdsges§ 


Heese 


sasttieerm 
sseenniee. 
Oaeet ee 


aere 
«. a9 
. 


Lede 


nani 


aaeee . 


cee peated 


in 
Hin 
Huntin 


3885533 TS 
pitti 


UiNmagEed 


Turbulent Film 


167 


The 3 Shoe Tilting Pad Bearing, 
c'/c = 1, No Pad Inertia 


Load between Pads 


Figure B-93 


ioe 


ee 


} 


rig 
aeedeesi we oeeee 


Aug 
sparsiesemansene 


TH 
Osea 


i 


8 Bassstttimasesatt 


a rel 
etesnss 
anieelit 


@e@ 
fa rit 


Tos 
CuUY 


ORO R SHIH! 
sotiiess 


UT 
s8sthumas 


=’ 


it aan Hy 


In EE 


OOGRINBE 
GORGES 


. 4 saver’ 
Hes a easeereis ee 


elect asieiiiian 
ma mamocenesntisiimecest 
ft 
SSBeranss UNHiNEat 


ait 
enaentn 


in ae a 
PT AURORA RAEI 


aut 
eegiitit 


Figure B-94 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionless Dynamic Coefficients 


169 


sancnai inten 7 oe Se 
A ne ng ees a 


=e SS re Sy ee ee ee -<-—s: 
EES Oe Ot OE ee EO ob ee Op Oe Ue ee hs ee SCCH SS Mss scesbane me? 
ee oho tks ae a Oo ones on aoe linseee tee iz i a 
F259 SEER SS Wem TCedd ES ee © OD oe ae Oe ieee aw aan ceded edge Be \ 3 
hbd4-debudpap-aht--detdidhdh toh dekh hee am 
Hdl ditpad oh Ame dehed Schot-d-4-4-edoed choad isebease: 
Bett ttt eed et 
iddtbgutAududed Adbb od 14-4 hod LLL AL LY rete hd 
ies tammemouiiseee diidtes GBeeil AB AUNTS 
By Ca ts Ww A 
BxHITISSS SS SSS e::: tat tig v } 
+t*steteseeowm a ae, tae iii Ks 
S280 OR tH faane rite d 4 
288eh oe mem a oe sen +a + 
Sisaeaoe seek 
ti-4 


. 
: 
* 


8609003) mn Oe leR 
99070:t1 RD marme 


Veta s ew tithes +4 
wigs Seawe bmw icele 5 
tpeSEBS* SS ices bi $44 

.*Sa 
s 

. ~ = 

.. 


amie 


Z: 
rs 


or Co 


a. 2-4 eres ee ee ee oe ee 


eS 


+8% 


: 
= 


Bahl 
aie 


= | © 


ee TA 
Ala 
at Fait i 


igs 


eee 
ine 


PT e tee eititit ttt teh 
Prete t Hei tt tied tet ts eed 


SEN 8iced Sekar es 


r] 
nit 
Ra 


List3gas 


maneauitt 


Dynamic Coefficients 


Dimensionless 


Figure B-95 The Plain Cylindrical Bearing, Gas Lubricated 


170 


ae | 


Bend 


esting OG 8e6s31h: 
0006s Sw OH OSSeErt: 


BO UULUIMNB OTe gett 
QRORRUINNBO REE 


LE Te 


GROTARINIG ROA GN Uii 
ST 


in prints ' 
mupene 
BeAlititaw heed Te Ln ALS aaa a 
ESLER hee eB Gnd bell ERE EE HIUI 


HCH etd d SRAARIBUOENUEH 


saiiitiesesmmeanscie 
99500)}) Gem mwa anand 
i eemene eases Beant 


" 
: 


Figure B-96 The Plain Cylindrical Bearing, Gas Lubricated 


Dimensionless Dynamic Coefficients 


171 


Hi 4 


as 
esses sitie 
eenges 


Heil 
tH 


ims. 
He 


Ssiivia 


Stsasttiiie 
gnansitive 


4 SEi2 
==s 


~~ 


ee es a a 


Nee 


SsSSSses itt cis : 
6 SSSSses Ss * attr Hts 
S=s==25 Ht sein Ses eessiiteseseinan 
ba dnad eae ¥, 
Pp) as an on oe Sen ee ew io 35:3: mma seessestaesonsscss 
meee am LJ aaa O6n98t2'h mmm e. Ht wenete Or Ronis: 
ttt e a8 me Stitt tte nena He aanesettn 
mw RA ee He these ae o8hei a Bane enehe: 
apt ee £09 me s16i1}) mimmeeeeastamumnenetennt 
Bees wensass var aeee aoe 
at _—4—_-—4- 2-4-3 4 ae 2 & 
SSSSSszs2i: HH 
Sapencsess: a 
BaBaszse::: 4 == 
ae a t+ amanae me 
2 a2 Gey Ge Ga ae setithitweee Hy 
Saemaaseaat a5 
maeease seeiti i HH ss seni ‘ss 


== mesasesit 
==ece 


FE Poe sth ea aan 
Z=582 Hi Eessr seine Sissi 
SSS sii SSsssssiimstini 
saazae siete soceriieenstee 
Smmsearsiiimiiiateiin sasstiime sails 
an =====: 5 dete ge cer 
ao ze = eae 
a i Suossesecess 


if imesviil 


=sS=SS== == ‘ 

= = (fa 

i= = A 

aa aewee : 
tiieese F H Nimes Hh ‘ 

aun j Himmmaness ili 

senas i HE Bee tetuniignntnniiR CUENTAWN ROOT 

TLLIT TT e A ETRE HHH eT COUT 

I 


2 & * * Sag 


Figure B-97 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionless Dynamic Coefficients 


172 


AO Nae a Yn Nate wees tee 


Sued ce es ee mee al 


Aight Sssless a SSS 

Sesehisisime aaeshidi hac ae cee 

tht oth iamee et 

rt ahaw ra a we at ae: 

Bacsesuiion r nagae jecess) ou am an emes 

Sescerucan iG { tt | HH ater aomoee oo oes 

ie geecttiiangia Tae" 6" iaganse win i-ttid 

SUUSIUINMMSUELEIHUMEE Ui cat te NideaegyEmeEnS 

3 BERS SSSR ssi SSSssiss SS SSSE 

zs Stim sssceuissessiittiuesne 
Lj 4 4 

3 Sesiiniudcssutivecsss 

aan a sousssssiawenaee 

men aguesorimaneneas ® 


st 
H 8 
see 
esses 


_— 
i — 


Et} SRS 
S32 


Ronse 
SS55533 
1] on 
sees srs somovaserscs: Ma 
osees os Saasmesces 


Geeae: Prepaid 
ad 


dd 

oe) mw Oa wees 088 
srvima we asanee 
HINMASR GI 
sumese ah 


ft at Ht 
wanes (awed ana 
onal | ttt 
MC 


Ptotiet ith 


ay 


i 
| 
Ao 
Poth t ' uy 
l gutlititii Hitneeenae PUTT ET) 
° 2 4 ie 2 4 a | 2 4 6 8 2 


10"! ae to 10 


S 
it} 
PS) 
~] 
Zz 

>30 


Figure 3-98 The Plain Cylindrical Bearing, Gas Lubricated 
Mueensionless Dynamic Coefficients 


173 


—~om 


i 
== = 4 ! 
== 3 mass 
Se cieiioaieaal oe tas tn t 
_— a om a eae - | 
aed mm Leenhancelpuinatcinal 11 aw i 
=a on itessnss siati mance rt rr 
oan se oe saaneess seme ioe ot ' 
mam 2a aoa sees Hert met He i 
men sa tiimesess tse eucer  dudegudld tt dg) 
=n =ain se Hesesiesinee asses: Hesse 
—4— 4 =. 2 re ae ec eses 5 p—4-4-$-$-¢ St tt waa eemente ‘ 
=a san SSSI ISS SSSSS a Sssssssss ; 
a ane weeen 3 ; Rccal + 4+T4 Siiimsecassres i 
CLS A i+ 8 2CPH OOS 0409 CER: waaanes a8 Meee nescat 7 
 ctandpmand ae se poe manne ea aawenneh ' 

si peaseesin === SRaseeess ae meneenel 

1 BBOsen: ht 6 Reseettstimase oenset 

sR eeeasisiiene mansnacesthioncenus 
ase sana MMAosreestt ma pe niet 
pie fr cree 

===. ttt tt ttt tet maednntily 

i=ms SSE SESS3S223 ii Sssssssete 

= = 4-1-$-44441-.4-4-5-4-44-42- 2h0 

= eA Da or i ee ee ‘lan natpapanteuchend dodecyl 

= Ee ONE ae ad te EE 8) MO ARO RE ee LEE dead b-ddnt dy-dpdeod 

nae at rn Sits 

‘aap GED cn eis mew Baese seri maneasnsen 
eam 0.00 oe Maweannsain 

Tt. 1 1 4 ewes wussseun Sees seen mame passin 
yy ene ee Bue wane npeehuop-depde det itu 
=S====5 SSS gHHEEES FERSSSs eee 
p—} —t—-} Ft +4 b+ —+—4 bt ——t— | ao seentenen wa arar 
an (or ae ae sabe eo t+F++ 4] were 
ES=SS==5 S== 33 H ssssctBni 
oo ir rot MRORS $beh a ae RENE OFe! 
a ane oan ieee Bees Pitt tt + ene on oe eaete) 
Hit aR a | maaeasesn pestsi 1s 84 81: 

Yt HH HHH Sndenutn thiee u 
sate HH at YY 


a 


nut me agatItinn zane , 
SRRRSASAHININIU SINNER Seta as tie 4 


il 
AE 
Beil 

a 
it + tf 
He 


=sSamsss St35i 25 SNES SS HIS3Hi: 
ses PR} eae pt tL 3: 
hoki] puget tt et 
+} BS SS SSB Soccs ien, =. ease ‘: 
A HOS tt sha 
ma wage Of Bums etBn,.:: 
tHE Ey FES | SEEREEH 
Sian ae ESE TEE 
se ietid) ines Sen sass eters 
seers Gis cet Beeee 200 ws ws Pease 
; | HE eugene shee ee w 
He uw Baralintineel wo 
im BORG rH aa i 

CO corr 

i=so8 Ss: 


Le}, iit 


eRe 
rin 
satil 


a je ain 
MBBS STON 010i 


a 
10"! .2 ee | ser 2 « 6 oe toe 


= 
e 
eo 
e 


Figure B-99 The Plaic Cylindrical Bearing, Gas Lubricated 
Dimensionlegse Dynamic Coefficients 


174 


eaaae SEPRAMIINNBEOEOLLII iia 


ons 
aie HH eageaitii 
RTE 


COT CA 
ol COC CT TT rT CECH TT 
io"! | 
re l@rpN 
Pa 


Sr ttt 
Hettrth seas 
enasnsem 
pastiisiase 
serenieen ime 


acest meme neg ss: 
UNNI OAR OTENNT 
RUMUIBROE GARNI 


Sitti nears 


ae 
SHIHIUUUNSHIUNUICOmBOREEEIINIOESIUU tin 


pF 
Et ttt tt tt 


Hitimean@meoesssie 


Pai 


nme 


Beer mea ai 
SSCS bse. 8 

ret te 

+4 

seteue 

eaesstim 


mOnasiiiiummanertitt 


on saauessey 


Figure B-100 The Plain Cylindrical Bearing, Gas Lubricated 


Dimensionless Dynamic Coefficients 


173 


0 Gees ree HSSSsineu SS Sesss 
J=aaesass pier HHSsssiusss seseseeces 
+--+} 5 haedeedbadh Leb pdbomdeadoabeddebedelcoh dbnowhmmmbedomcdsArchudedechd Fh WOME eee 
—o—o = Ht Pt ttttte ttt 0 TD Saitama maBenenet 
3 SseH ESS RSUSSls) meee ene Bamesscer craae measesy: 
iaseinnanmmaesereimseenini: 
IMGSSGECIii Guan mameasesans ea} MOO Senses 
st SS Sssestis! ssa asesseises et wee we wr anasa: 
PBS SS S555 SAS Soss sess Ss Sesisthi 
ed ee Sr ee T adnelieadhandpedpetiyd.cn.dehdehuradiaaleraduaiend Aydt 
htt Steed hdd Ghana mawensee | Lo divdbandbedaMad hae 
WO OS 60822111 am eee eesate asttmm rite 
4 5254052 Gn aan eaeeeeecee 2e aaseees 
ot a a Geseer 
i] peat a a jee aaqat 
Secueresit ma maceissinimanecie 
SORE Coo 
BARRA Pr TT 
mages mss ee STM RO TIN 
3S =S= ; Seses: 
8 SSSessseese : 
om b+ —+ 3 L, ee 
if ne tt z3 i ttrt mae neeseisie 
5 ar 2 ons OA OR BES 8, . Sc a 
Seaas Baer aitmisesetinaitan. “wassnnnisaseii 
ro ind - 
juname iapenr.* afouinriiposscinsitesin, ‘“eesiltuimeeseatit 
= = f j Seceist pS 
== PF FF HE: Bassese "areses 
Sica ai fist 4 Ti. Lt Oe “@eeea 
Raase ? 7 aBvi8 RO Os ers) ame ean: suaeee. Sedat 
-_ a eae B6dlels Ra mwRe AAS eeens . Cae 
So ay oe f Pee O02 mmMmane i j+t+ “s 
toon ons ous ons oo ooo Goenenmasccosrin ettnsemenesstinitmsas 
anene ae meemmaenses noe Coo c ert 
4 aaa Mmeeausenes nsesl smesnstsrnuinaatttii 
EEC Sesiiiiiinal ust meaattaataatia 
10 seeeere i= B255551 RSS Si RSS SESS Ist S35 555535 
e —— 
i ee Eee 
}—+ ttf +4 
7a Gee 
sae= 
then bn fp —pneed 


2338s: SSS, 
hth tet —4 


Figure B-101 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionless vynamic Coefficients 


176 


SUT BO ONL 
UT 
TO TT 
OTT 
Attiim@esansinn 
ptmipeaakes 


ae HiT ae 
woenese = ere? awe |. 7” 
sss Sea Heduee ie remndeassisiimecses 
weckesee eer t ttt tt tte 
1m ne Pritt ttt ttt meme nentnt: 
a60se: 1. peace se me anaeninn 
beet wr ences: 
annsst cones A004(05iM0000I0IIIi 


rT Ci NLL ULBRUEIEN 


detrei: 


Saas eee 


= it 
= i 
as Ht 
a tt 

a Bueeriil at 
a QOUULIUN MR aREth 
ot LEC TTT 


Figure B-102 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionless Dynamic Coefficients 


177 


hus mssasssiiniseais aon ease 
We a8 61060) tee ett ret ty oan ee. anes 
sesieeuisesasisisisianmmmesasstiumesssssins 
O's Boaesesss!sian ance suensccsiscamsa Lesa boed 
a 8 ‘ima actttitil manmen susuiitutmeaessatt 
eereseen a: S00iimanas eanse we anne 
i 


s ea 
att tt ese ence: SSSR ESSe 
Seeessgeestiion 
ona 


Beoos 


im 


ae ee iy 
sntieetet oe edt ans oni ROO OE 
sities since ‘ 

ype 


af 


, 
ee 
fe 
a 


hi 
if 
E 


; 
UUmmees 


ae 


id 


am avarini 
riuMl dttittion ti tat 


tar 


| 
fF 
aH E 
y 
BF £328 


isiaeasii Letra 


+ 
‘ 
d- 


tt 
BITE ct t 


: ad Ht SH 

Pt wees HHH Ra Ba aN 

Seadaniy ail aunt jes Hit ee 

0 OE CROHOECUUTIIIGSQUIIIII MBGOCOTITIIIGGEE Hee se Tt 
° 0”! : 4 6 8 ; 2 + 668 0 2 "102 


l2 wu N 
Ae ome re 
Py 
Figure B-103 The Plain Cylindrical Bearing, Gas Lubricated 
Diwensionless Dynamic Coefficients 
178 


aagae ee a 


serene 


seapans HigGuacegtnttananes 
Huageniitioas 
Onis) a 


uty ie 
eel ea OT HH 
atte each TORE TT 


Q 4 68 2 ¢ ¢ 38 a 4 6686 


Figure B-104 The Plain Cylindrical Bearing, Gas Lubricatad 
Dimensionless Dynamic Coefficients 


179 


: Serer er rere 
i ama aessasqesienesants 


SS oH sein: HES ses 
tt t+ tH Sings Sttsstieese: 
sie sn ns Sus om auton ene sgeset me buaessst\msomme 
=o aaneaeanea imasesisiii me 
HIS tt 
mescsssinmecas time 


ocnnaas r 
oe BB Seasenses 


Seessiteisie HH 


am 4 HH swans 
EH tt 
ees eecrsinis mannered 
erreimeas moaesit 
esceti MGABarnenD 


BSacss i) HH Hit eae 
ee 
BESSRESIA CEI pate tee 


tt aan 
Sasesas siiuiniesss 


Sassssiieeiauk 


SESS: VS SV EBESE? 
2: 33 3 


« “a. 


S23: Ses 
Ssfi2: tee SSwes... 


Hf 
H 


tt 


tits 
}. 
il wanet 
a 


ttt bos 
rooe corel 


Ho 


Si SES Baees 


t+ ee te 
ae t. 


[: 

4 
BL ‘peur: 
Bi phe 

ever: 


Ss38E5 5; 33 i @aesiisnarstiinn 


Q! PH al 


2 4 ee a a oe 2 4 


Io! | 10 he? t ! 
| 


Figure B-105 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensioniess Dynamic Coefficients 


180 


sa aeangmaaneiiinat tae mom cece ee ee 0 tr errno magnate + ooerarnente 


= = 


[ 
{ 
I 
z 
i 


ae 


SsssissseSSTS seis: 
ms SSs552523:: scram pomeckises: 
| Sede S20 a8: 
HSassistiineces sists ih 
siatiisa 


asus =sss Sacareseumasgsegsaie 
aeastiiimemaas Bessesibimn adnan Ha) 
001 OO 0096111 Sw anes OO Se skeet ms we Rs 88 8st) 
& SSGRsttHitimaassteitit 


sine Ee 4 
& mnae ORCGHEAIHIMUGUOIHUIIU 
| aR HALEN fl 


Pritt Hoth 

ae araeeniiien 
SCR R ae cael) 

TAO BUOAEIUi one 
IHU ROBB ANIL 


INwaee nes) 
Tt med pace 
iat tnte 


man a a we 
BS333; Poaet aes wee wa a a. 

oage am onan re Dees 

SSSs3Ei5S552 


omnes Ess + eee 


ostaim 4 
CO 
Begesiiiias Hy 


8 
LU 
t 
t 
1 
! 


an aan au a8 HTH REEL RRL 
ot Oe CHUOGGUTAIIUIGUQUGTIIIIGUMMOSGKIIIITIG 
10°! 2 4 5 8 2 4 6 10 102 
re l2 mes ay 
(i 


Zigure B-106 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionless Dynamic Coefficients 


181 


eee ———— ass sneer 
— S SSSSS SST SESS Hi = 335 S==re Sotcict ease seatcss 
a eae ttt ttt ttttttt rt + +t 7 are ae Zs Bu a. PF -4- aennes TOS RR REGO RE Es 
Laman pp 2 eens s t+ rs CD CE SD Gib =P =P ee esss: ms oaese => a oF ORS OS Res OF a te | edd --debil 
SSSSSSe5=55 HH Saas See eeari SS So 505i2it SSS Se ssssesineossststse: 
8 ttt ttt oe end BSwasccessei aca Saasasets: aden ceana ea, YL Lb bdddeod-old OOH 85.09 
+ tt PP oe Ht +4 Sceasstisiase enanar Hy te tte H+ 4 @oasann 
oe Mmm Sen e Beaganesei lad sanenenses: eens eo mmemnansece it OS Bees het: 
tot ttt eengae HHt-tt +44 eeeess UNSOSSESSSSImeammensresss ti maatenn 
_ Suastsitionm SHAKTI MMAR TET MeN PRaAAsAte MAE Ae ATETI 
pt [ttf —*-* $F $4 3353 tit -t-F-$- 4+ ttt —-$--b—-$--F-F-F-J ake SS255 5325 
= + t+ ++ a. 33. sce. eS SBSes80 Pe Fh P 4 na Samsaaern 
a HE St = Oh cae ETS I aE Sosis esemMeasesn capes wee ee See ee seesae a 
EE: EE GANG CRG OES ME One VER ED ER am oe WE eae Se rit ttt tty 01 ae a Oe Ek ae aS O35 DP tee Een etat 
Ry ae RD BE ey Be Gh Oe oa ‘Ree Fee Sore de A ee ee Benan eH HE 
~~] =a aneeees apeem y= HEH tt Pi To 
am aD OE Me an a aes @e CBR eues ii 20 GP 0S GSHSHS 6001 eM MA Hs Hd8d81 
tm Gre te Bm ms ee T] @pe*ey ar iBSSeasessitinnaratmemensststiimeeastartt 
p] MASS ee ie it Ht =e SOSCAl 20) Ma SS SSeeststimene oeuenie 
Pt HH aes Sasisititiiewstmasaaitit Nmomopentet 
= oonnoeeent Cott tty Manaaccins mee SASSAsei BERR RAALLL a 
Ht SRST HH IIUSSELRLIIOnaEs Hoth HHH HCH BCH tt 
ne mmmOOnaesti HiME EAL SET ttt SH MOB 

ane ee =s VSSisssSSSSsss ss ese 

e SSSSSEF223 t+ tty Sates si es omnes ches! ees am 

co +t $+ oo tt aa 


Be esse sista 
SHnsSesasssrssitusssrsiiiii: 


—_ [J 
22 a Oe On oe a WMORa ne ntsii 
Pia inencaas meweNLE imanpad 
ZRSSS5 = 
EG SE CSS: ONCE ee ae 
Sas=5 Se eeesetcit 
PY ce ov eesti GAM Bees sasie 
RE OR es Piet tility 
Smee m [ . Pd Ge 
S2eaaee Ls Po TT 
a oe eta Tee ed 
CCT U1 Do, SHE MT 
= = 4 BS BSSSS Sasss ese SS STss HH 
ry 
eet 
6 


Hey 


ns Ea os 
OERLUATN MARDI UIEI 


ry ise 

ge HituhEs 
‘ =e es = i 
seeesiiiieiiiiiee Serie 
3 HH HCH Hiiimanes aeee eeniim 
10 CHT tr wu PCO 
SMOARRAtISIIINMUGNOKI Ce i ACTH 
TTT TTT ie 
ae oaaarnait HIM OPAUHLH mani He 


Figure Be107 The Plain Cylindrical Bearing, Gas Lubricated 
Dimensionlesa Dynamic Coefficients 


_ 182 | 


p++ HEE setretiy 
1 aoe aes smansenssrt soounasst 
itt Se ee aaisi 
OR ASSES) rH 84 oem: Ha 
Ha earsse Sime SEnBtstey 
pid SOOVTIII MAw|ss sit 
EHHiL eeRmUE aun AniBaseRtatin 


- 
se 
me 


eesas 
i 
* 
e 
? 


= 


a 

Leda 
+ ‘ 

= 

= 


ret 
er 
H 
ae 


WOE OS ESEEE 
MmResus , 


Sen Ma ee | 


aaa ae | 
Tag 


Serna 
Socetiiueeccss 


Sqssiymes: ep pestis 


ms meastaiit Po auininsas 


| 


a 

Tt 
ime 

vES 
HIEES 

iim 
imme 
Hees 
saa 
‘=== r 


= 
Ld Se oe a 
+ 


i 


aoasese 
ai 
sons ST 
Co 


ssesasatiia ry 
MTT imaaan Saeessoneian Men ernia 
fA ‘ asst tH 
t Ht a Pa St 
aun HN 7 AiR Se mane 
ND. @r/a7 Tal ea 


Te ee a 


iat 


ce d 
by as : 
saber 

, MOan arenes, i 

Li buded- pb LE tit ' 

; s8SR SS SS erty 
ORS Oe Hh OO M8 Od Rss, 4 
SONF em DS ePOL Ty) 
O06 AO OO BASS Badr ‘ 
ispsssssigie 
esi ma ens: tts 1 
ateim@e a Ht 

eeriaaG E) Ht ; 

' 


pera rc ae 


= 
= 
= His 
e 


Hi 
BOSUAtifnimmasssertit 
COM 
(iat 
QLGCEITH MAGA AHINH 
RTH Hd Hitt HORGUSN MORAN! 
MEMO STRTT HMA H RANA it ut it PERU ROTA it 

oh PE CECLCCETETE En CCC ECTUUIM TR URN ALSt! 


e ¢ 866 66 Q 4 ¢6 2 * * “0 


io”! a 10 
. G | 1 


Figure B-108 The Plain Cylindrical Bearing, Gaa Lubricated 
Dimensionless Dynamic Coefficiente 


PEC l Cth a 


183 


Petty ett Pt tt ee ee eee 


SEtssilihe SSsssesenes 
= 


ne 


Hit +4 
BE 
Hise 


Bevesnriieas 
ae 
eam 

titi 
NHI 
pees) 
Subsea 
3 
ty 
Het | 
u 
Reid 


Lyi Ht 4 aoee 

nant vameeas 
meu. eer 
aoteee 


f = 
BUSdSs. were c—_ 
Thee 
a oa SSanaseasat tt 
= ssasss =a 


ae? 
- 


333 petty tt itt hss 
S22: Ewes! seesess se 


Dated 
— 


He 
bt 

HE 

an 

IT] 

: 


rnhis 
mmeanesss saisinn tH ‘iwas 
meneeens as 


See - 
t=a 
338 
22 

2s 
22 
f 
cine 
a8 
= 
fi 
+44 
Ses 


. 
Lee 4 


oor ape Deen sens, 


Sesiiit 


P,LO SS asledes sceceeisttl a : meee rf [ | 


ee encccsnian 
ToT 
iy 


10% i it ete Ett te Ht 
Rana Hott a EHH | 
BEC TTT : 
(EEE ae he 
e Tea 2 4 6 86 2 4 6668 10 2 4 @ 102 | 


Figure B-109 The Piair Cylindrical Bearing, Gas Lubricated 
Dimensionless Dynamic Coefficients 


184 


CE Ree aks 
cw See 


re ne : : Oko Baia ‘ 
cme oe LEN TER eT wee etl RED 


d. 


Gas Lubricant, Hydrostatic Bearing 


185 


UoTESTmpy suRTd aTBuys ‘o> = ; ‘8a 


ae | : a 
Trt tee ---- - - I. OS 
Nich : : i ; 
{ thes 4 


SseuzsJIS sseTUOTeUIMITG 
FAveg SPD YFIBISOIPAR aq O.{-g aan8q_Z 


= $020/7 — NOISSINGY INV Td- TONS 


: : } poe 
Be tp eee ne oe tee 4 ys fe 


160295 


y ¢€ 


TOnNo mn bal 


d) 28st 


an ‘d- 


28+! 


» Be) 


186 


PS Seer Gee ae 


. 
rt 


5} 

+ 
4 
Sia 
7 

{ 
+. 

4 

4. 

+ 


ee ety 
Lee & Pe ee 


+ 


pasta OP ane tn ep me Beg ee DFAVIOIPAR —4E TTI-fl BINT yg 


oe eee one 


. so 4 aS eo wane awe 


ho oe 


ea 


eee 


= @ 


187 


: ~ wt aie 


1 2 

i 4 See ig A 

\ We ey EG 4 
panne ee HH 

i ew 1 ytd 

Satan tt et 

‘ ee can 


| 


tht ppteote 


' 
a 


aay 5 
a ao 


Date 
‘ 
ot a a pas oie 
ty 
tt 


188 


OMISSION: - 


SINGLE PLANE A 


Og, ke -3, Single Plame Adsiesion 


Dimensionless Damping 


Figure B-112) The Hydrostatic Cas Beari 


‘a SSIVJJI25 SFT TUCO}; sUSMTG 
| uoysstepy suvig a[Zuqs ‘{, = ¥ *Suyzssg sey 9798Q801PAH BUL CII-q JaNn3Taz 


00! 01 


| 
‘oO 


eee eae a abe seated 5 aes ’ t TS ge os t +. > = a'¢01 
1 ‘ ; | i ! 4 t i 
ic cal i i 
} - -j+ = 7 I + 77 z 
' ' t ! | ; i. ; 
fy oe ae ease (" 
“a1 7 _Noissinay INV'Id JTONIS 2 
a Nidan Cae ee paeiees | + 
lo 
jo, = 8 
+ |— 
wire] + 
z oo 
baa ~ 
—-—{¢ ‘ 
- : a) 
pce pecs ofS ’ 
ears | Phd Se 
6 - 
; ‘ oO 


t 
arte aeranr oc rs] 


Meme a ie a eee eR | 
6 2 wegz9S h € 2 t 


teucs & vr & z tee. S vr € z 16e¢ 


ve ; Peres 
fre ae eee ts 
roe eee - oe pay 
i + : 
‘ + + ¢ > 
Beem Be eth 
& qbeeneesteune 
senves 
2e8 He 
e © 


pee.” a aE we # 


Ae ee ere BREET PEASE Sess? 
ES OES 
Rei aan base EEE aioe y a 


190 


Dimensionless Dampiog 
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Figure B-118 The Hydrostatic Gas Bearing, 
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Figare B-122 The Hydrostatic Gas Bearing, 
Dimensionless Stiffness 
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OF A_ROTOR-BEA 


Discussion 


. 


A vctor supported in fluid film bearings is susceptible to an irstability known 
as fractional frequency whirl ("oil whip"). The instability mechanism derives 
from the hydrodynamic action in the bearing film where, at a certain disturbance 
frequency value, the film looses its damping ability. This instability frequency 
is an inherent property of the film and depends on the operating conditions of 
the bearing. It is normally expressed by its ratio with respect to the rotational 
speed. Thus, if the frequency is denoted as vy radiana/sec and the angular rotor 
speed is w radians/sec the frequency ratio is v/w which is dimensionless. The 
value of the ratio is a function of the eccentricity ratio € (hence, also of the 
Sommerfeld Number S). For high eccentricity ratios most liqui@ lubricated bear- 
ings are stable such that = is zero or non-existent for approximately € > .8 As 
€ decreases ¥ increases Satdiy and levels off with a value of .5 for all € < .6 
to .?. seacas hydrodynamic instability is also called fractional frequency 

whirl or "half-frequency whirl". 


For disturbance frequencies less than the instability value the effective fluid 
film damping is negative becoming positive for frequencies higher than the in- 
stability frequency. If the rotor-bearing system is disturbed away from its 
equilibrium, it will seek to whirl with its lowest resonant frequency. When the 
resonant frequency is higher than the bearings inherent instability frequency the 
disturbance will die out and the system is stable. On the other hend, if the 
resonant frequency is less than the instability frequency the negative film damp- 
ing will cause the initial disturbance to grow and the system is unstable. There- 
fore, the threshold between stability and instability is determined by the re- 


quirement that the instability frequency equals the system resonant frequency. 


The resonant frequency of the rotor-bearing system depends on the flexibility of 
the rotor, the stiffness of the bearings and the mass of the rotor. Thus, to 
define the threshold of instability by means of genecal design charts requires 
that both the rotor and bearing parameters are included. In the present design 


charts, the rotor is assumed to be rigid so that the rotor mass is the only rotor 
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ice ee ee 


yee eee te een 


parameter. The additional effect of rotor flexibility and support riexibility 


can be estimated by calculations as demonstrated later. 


se of the Design arts 


The design data for the instability threshold are given in Pigs. C-1 to C-17 


covering a wide selection of bearing types. 


a, Liquid Lubricant, Laminar Film 
Figs. C-l to C-7 apply to bearing types with an incompressible lubricant and a 
laminar film. As a function of the Sommerfeld number § (defined cn page 27 ) 


the curves give the dimensionless critical rotor mass: 


Dimensionless Critical Rotor Mass «= a 
upi(®) > 
where: 

D - Journal diameter, inch 
R - 7 D, Journal radius, inch 
L - Bearing length, inch 
C - Radisl clearance, inch ‘ 
Hh °° Lubricant viscosity, lbs.sec/in* 
W - Bearing load, lbs. a | 
M - Rotor mass per bearing, lbs.sec/in 


The rotor mass per bearing may be computed as: 


2 
‘: -[- istance from bearing to of rocan] - (Total rotor mass P88 : 


(Rotor span betweer. bearings) 


Thus, if the rotor is hcrizontal and gravity is the only radial load, then: 


Me E (horizontal rotor, gravity load) 


. where: 2 


g = 386 in/sec 


The procedure for determining the threshold speed is then: calculate the di- 


mensionless rotor mass according to the equation given above. Enter the appro- 


priate design chart and read off the corresponding value of the Sommerfeld 
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) 


number. The speed corresponding io this Sommerfeld number. is the threshold speed, 


i.e. the rotor speed at onset of instability. If the operating speed is higher 
than the threshold speed the design must be revised. 


The charts assume the rotor and the pedestals to be rigid. However, the threshold 


speed with a flexible rotor and flexible pedestals can he calculated on the basis 


of the given design charts, Let the stiffness of the bearing support be K, ies 


and the stiffness of the rotor K.. te, The letter stiffness can be taken as 
the lateral stiffness of the rotor by Ltself when a force is applied at the CG. 


Make the two stiffness values dimensionless: 


ou cK. = CK, 
K,. - W K, > W 


and introduce the parameter: 


ite K.° K. 
K+ K, 


Note, that if the rotor is rigid: 


k= K, (rigid rotor, flexible pedestal) 
whereas for a flexible rotor but rigid pedestals: 


k = K (flexible rotor, rigid pedestals) 


Next, choose a value S of the Sommerfeld number and find the corresponding value 
of the dimensionless rotor mass from the design charts. Denote this value as 


M,- Define a dimensionless effective bearing stiffness by: 


—~ y2 
kK, > (xsi, ) 


Then the dimensionless critical rotor mass for the flexible rotor and the flexible 


pedestal, corresponding to the selected Sommerfeld number value, becomes: 


= f k 
= M — 

Z| ° \ k+K, 
HDLC) lex.rotor B 


flex.pedest. 
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Repeating the procedure for several values of the Soumerfeld number defines « [' 
new curve for the dimensionless critical rotor mass. Entering this curve with 
the originally calculated value of the dimensicniess rotor mass deteraines the { 


Sommerfeld number and, hence, the rotor speed at the threshold of instability. 


No instability design charts are given for the tilting pad beariag. All the de- 

eign data for the dynamic coefficients of the tilting pad bearing assume that the [| 
pivoted shoes have no mass inertia and under that condition the tilting pad bear- | 
ing is inherently stable. In practice, of course, the shoes will have some [ 
inertia but as long as the phase angle between the shoe motion and the journal ‘ 
motion is less than 90 degrees the bdaring will be stable. The condition under [ 
which the phase angle becomes equal fo 90 degrees can be determined from Figs. 
B-37 to B-48. Here the dimensionless critical shoe inertia, labeled "Critical [ 
Mass" on the curves, is plotted against the Sommerfeld number. It is expressed 


by: CWH nae [| 

Dimensionless Critical Shoe Inertia = y—. 212 4 
[pons 7] 
where: ‘ [| 
2 2 
Merit = I/R°, 1lbs.sec’/in | 
I - Mass moment of inertia of a shoe around «a jongitudinal LU 
axis passing through the pivet, lbs.in.sec 

R - Journal xadiual inch iy 


1 
The other symbols are defined on Page 27. 


Hence, in checking the stability of the tilting pad bearing calculate the 
dimensionless shoe inertia, enter the design charts and read off the corres- 
ponding value of the Sommerfeld number. This value should be greater than 


the operating Sommerfeld number. \ 


b. Liquid Lubricant, Turbulent Film 


Figs. C-8 to C-14 apply to the plain cylindrical and the 100 degree partial 
bearing with an incompressible lubricant and a turbulent film. The curves 


| 
| 
{ 
\ | | F 
l 


give the dimensionless critical rotor mass as a function of the Soumerfeld 
number S for a wide range of Reynolds numbers. The procedure for using the 
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charts is completely analogous to the procedure described above for the laminar 


film bearings with the exception that the dependency of the Reynolds number 


on the rotor speed must be included. The Reynolds number is defined by: 


Re = 22NDC 
y 
(for symbols, see Page 27). N is the speed in RPS ind is also a factor in 


the Sommerfeld number: 
_ HNDL /R 2 
5 


Hence, in using Figs. C-8 to C-14 compute first the dimensionless rotor mass 
(see Page 202), enter one of the chares and read off the corresponding Sommer- 
feld number. Based on this value calculate the -otor speed and the Reynolds 
number. If the obtained value of the Reynolds number differs signi fican:ly 
from the Reynolds number for which the chart is valid repeat the procedure until 
both the Sommer fd number and the Reynolds number yield the same rotor speed. 


The resulting value of the rotor speed is the speed on the threshold of in- 


. Stability. 


¢. Gas Lubricant, Hydrodynamic Bearing 


For the gas lubricated bearing the design parameters are: 


Compressibility number: A= Deel &? 
a 


Dimensionless bearing load: —a_ 
P Lb 
we 
Dimensionless critical rotor mass = ae 5 
He Xa) 


where: 
N - Rotor speeji, RPS 


Pa - Awbient ys essure, psia 
t : 


and the other symbols are defined on Page 27. The rotor mass per bearing, M, 


is determined as out .ined on Page 202. 
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Figs. C-15 to C-17 apply to the stability of the plain cylindrical gas bearing 
for three values of the length-to-diameter ratio. The rotor and the bearing 
pedestals are assumed to be rigid. To use the charts compute the dimensionless 
bearing load and the dimensionless rotor mass, enter the appropriate figure and 
read off the value of the compressibility number A. From this value calculate N 
which is then the rotor speed at the threshold of instability. The rotor-bearing 
should not be designed to oferate beyond the threshold speed. 


If the rotor and the bearing pedestals are flexible the threshold speed is lowered. 
Let the pedestal stiffness be K, ibs and the rotor stiffness x. ibs as discussed 


in 
on Page 203 Make the two values dimesionless: 
_ cK, - CK, 
Ke" Pop kK," PLD 
a a 


For the given dimensionless bearing load choose a value of A and read off the 
corresponding value of the dimensionless critical rotor mass from the design chart. 
Call this value M- Define a dimensionless effective bearing stiffness as: 


2 
vaN -_ 
Ky" 288 ™ 


Then the dimensionless critical rotor mass for the flexible rotor with the flexi- 


ble pedestal can be calculated from: 


| $s] oa a 
2, ,R, 5 k +k ° 
: Xe) flex. rotor - 

flex. pedest. 


By repeating the procedure for several A-values a new design curve is obtained. 
Entering this curve with the originally calculated value for the dimensionless 


rotor mass determines the value of the compressibility number and, hence, the 
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totor speed at the threshold of instability. 


If the desired operating speed is above the calculated threshold speed, it is 
necessary either to revise the design, or if this is insufficient, to go to 
another bearing type such as the tilting pad bearing or the hydrostatic gas 
bearing. The latter bearing type is considered below. For the tilting pad 
bearing no design data are available for studying its stability. However, as 
long as the mass inertia of the shces is small the tilting pad bearing is 
inherently stable. An estimate of the critical value for the shoe mass inertii 
can be obtained from the design curves for an incompressible lubricant given 
in Figs. B-37 to B-48 when the proper bearing dimensions and gas lubricant 


properties are used. 


d. Gas Lubricant. Hydrostatic Bearing 


When the speed of a rotor supportedin hydrostatic gas bearings becomes suffici 
high the hydrodynamic lubrication action in the bearings starts to be signi fics 
A hydrostatic bearing operating under these conditions is called the hybrid ga: 
bearing. Due to the hydrodynamic action the hybrid bearing is susceptible to — 
fractional frequency whirl ("hybrid instability"). Only limited design data ai 
available on hybrid instability (see Ref.19) and because of their incompletenet 
they are not included in the present report. However, it is possible to give - 


simple rule for obtaining an approximate estimate of the threshold speed, name. 
Threshold Speed = (1.6 to 2) + (ist Critical Speed with purely Hydrostatic Bea) 


In calculating the lst critical speed the hydrostatic bearing stirfness obtains 
from Figs. B-110 to B-124 should be used. Since the rotor for a gas bearing aj 
cation is almost always rigid the critical speed may be computed from the 


‘equations given on Page 235. 


Besides hybrid instability the hydrostatic bearing may experience arother form 

of instability called pneumatic hammer. This instability is in mst cases in- 

dependent of the rotor speed and is, therefore, basically associated with the 

design of the bearing. It derives its energy from the supply pressure and is 

analogous to the action in pneumatic tools. A design of a hydrostatic gas bea 

ing can be checked for pneumatic hammer by means of the design charts for the 
207 


dimensionless damping, Figs. B-11l to B-124. Whenever the damping te negative 
the bearing is unstable. Hence, to use the charts for the purpose of checking 
stability calculate the restrictor coefficient A, the supply pressure ratio 
Pi/P, (see Page 69) and enter the design charts to determine if the damping is 
negative. If this is the case the design must be modified or the bearing may be 
made inherently compensated (i.e. the feeder holes are not provided with orifice 


reatrictors) . 
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Figure C-1 The Plain Cylindrical Bearing, Laminar File 
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Figure C-2 The 4-Axial Groove Bearing, Laminar Film 
Critical ged Pca at Onset of Instability 
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Figure C-4 The 50 Degree Partial Bearing, Centrally Loaded, Laminar Film 
Critical Rotor Mass at Onset of Instability 
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Figure C-5 The 60 Degree Partial Bearing, Cercrally Loaded, Leminer File 
Critical Rotor Mass at Oaset of Instability 
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Figure C-6 The 80 Degree Partial Bearing, Centrally Loaded, Laminar File 
Critical Rotor Mass et Onset of Instability 
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Figure C-7 The 100 Degree Partial Searing, Centrally Loaded, Laminar File 
Critical Rotor Mass at Onset of Instability oar 
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Figure C-9 The Plain Cylindrical and the 100 Degree Partial Bearing, 
Turbulent Fila 
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Figure C-10 The Plain Cylindrical and the 100 Degree Partial Bearing, 
Turbulent File 
Critical Rotor Hass at Onset of Instability 
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Figure C-1l1 The Plain Cylindrical and the 100 Degree Partial Bearing, 
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Figure C-12 The Plain Cylindrical and the 100 Degree Partial Searing, 
Turbulent Fils 
Critical Rotor Mass at Onset of Instability 
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Figure C-14 The Flain Cylindrical and the 100 Degree Partial Bearing, 
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VIT 
THE CRITICAL SPEPIS AND RESPONSE OF A ROTOR-BEARING SYSTEM Ly 


Discussicn fi 
The rotor-beering system is a dynamic system which posseases toth inertia, , 
fiexibility and damping. The rotor itself is an elastic beam of variable cross- sf 
section and with a variable mass distribution. The fluid film beerings are rep- () 
resented by a set of 8 lirear coefficients: spring and damping coefficients, and {} 


the pedestals may also be flexible and have inertia. In this way the overall 


system may be represented as shown in Fig.12. A more detailed representation [ 
u 


of the bearing coeffictents is shown in Fig.13. 


The very complex system made up of the rotor, its bearings and its supporting 


structure possesses a large number of resonant frequencies. If any of these 


resonances are excited in the operating speed range large amplitudes may build | 
up and prevent further operation of the machine. This is clearly illustrated [ 
in Fig. 2 where tie amplitude of the journal in the bearings becomes very large 
upon encountering a resonance in the pedestal structure at approximately | 
24,000 RPM. Not only is this amplitude too large from a safety point of view 

| 


but the resonant condition also makes it much more difficult to balance the unit. 


Even though, as illustrated by the above example, other resonances may be im- { 
portant the most pronounced resonances are in general the critical speeds of l 
the rotor-bearing system. They are excited by the unbalance in the rotor | 
and the amplitude is only limited by the amount of damping present : 
in the bearings. It is, however, evident that if the bearings have too much f 3 
damping they sill behave as if they were rigid and,thus, would be unable to | 
dissipate any energy from the rotor motion. Consequently, the rotor amplitude | 

| 


would be very large. 


To calculate the critical speeds and the unbalance amplitude response for an 

arbitrary rotor configuration supported in fluid film bearingswith 8 dynamic | i 
coefficients requires a computer program. In such a program the rotor itself 
is approximated by a number of discrete masses connected by flexible, weightless 
shaft sections, see Fig. 12. The shaft sections have the same stiffness as in 

the actual rotor and the actual mass distribution is appzoximated by the discrete 


The idealized rotor can be brought to represent the actual one with 


: 


masses. 
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as close an accuracy as desired by choosing a arge num * meses stations. Th 
values of the bearing coefficients to be used «= the <a: -Lons can be obtaing 
from the design charta given in this report, see Fig. Bel to B-124. 


Although accurate calculations require a computer program it ie atill desirable 
to have information available from which the critical ape ia and the unbalance 

vesponse can be estimated. Such information is especial? y valuable when compar{ 
different bearing types and their influence on the performance of the rotor-beat 


ing system. The design charts in Figs. D-1 to D-40 have been prepared for this 


' 


purpose. 
Uge of the Design Charts 


The design charts are based on a simple rotor model, namely a rotor consisting 
of two heavy discs on a flexible shaft and supported in two similar fluid film 
bearings. The rotor is symmetric about its midspan and the unbalance is intro- 


duced by letting the two masses be eccentric by a distance 8 inch: 


The chosen zotor model can also be used to represent a rotor with a single centz 


mass. Hence, the results are equally valid for both rotor configurations. 


The rotor stiffness is expressed by means of influence coefficients. Thus, ay 
is the rotor deflection at the location of the first disc due to a force of 1 lk 
at the same location. Similarly, Aho is the derlection at disc no. 1 caused by 
1 lbs force at the location of disc no. 2. Due to the imposed symmetry of the 


Introduce the 


rotor the influence coefficient Gy> at disc no.2 equais Qy° 
definitions: 
a = {Su ~~ > (1st critical speed) 
; M17 % (2nd critical speed) 
1 (let critical speed) 
“te see Fig. above (2nd critical Speed) 
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The dimensionlese design parameters used in the charts are: 


Soumerfeld number: S = BND (By ? 


a 
Rotor flexibility parameter: pop # 2 
Cc 
Speed ratio: @ 
bac) 


Normalized rotor amplitude: - 


where: 
Journal diameter, inch 


i D, journal radius, inch 

Bearing length, inch 

Radial clearance, inch 

Rotor speed, RPS 

Mass of one disc lbs,eac/in 

Bearing reaction, lbs. 

Max. whirl amplitude at location of the discs, inch 
Eccentricity of CG of the discs, inch 

2"N, angular rotor speed, radians/sec 

Critical angular speed of rotor in rigid bearings, radians/sec 


3& €E ox BF Baerewes 


The rigid bearing critical speeds can be expressed as: 


1 
Oo = 


n Mc ‘ 


which gives both the first and the second critical speed (rigid bearings) depend 
on the value used for @. The unbalance induced by giving the discs the accent 
ricity 8 is such that for the first critical speed the discs are displaced in th 


. game direction (static unbalance) and for the second critical speed the dis~ 


placements are in opposite directions (dynamic unbalance). 


To apply the design charts to an actual rotor-bearing system it is necessary to 
know the values of the rotor flexibility parameter p and the rigid bearing criti 


speeds @ However, an accurate calculation of Os especially for the second 
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critical speed, is not readily accomplished for an arbitrary rotor geometry with 
eccess to a computer program. On the other hand, if a computer program is avat! 
anyway the design charts would be of iittle value in investigating a particular 
rotor-bearing system. Instead, the purpose of the design charts is to make it . 
possible to compare various bearing types with regard to their influence on the 
critical speeds and the unbalance response. Hence, for using the design charts 
it is not necessary to have accurate values for the rotor flexibility parameter 
and the critical speeds of the rigidly supported rotor. On this basis the resul 


obtained from the charts will give relative values. 
Tc investigate a particular rotor its total mass and mass mowent of inertia ard 


a transverse axis through the CG should be computed. Denote the total mass as ' 


and the mass moment of inertia as I. Set: 
E - ‘hee 
ul 


2 - Rotor span between bearings, inch 


where: 


Evaluate the rotor influence coefficients %,, and Gio ve respect to position 
on either side of the CG, i.e. position 1 at a distance 3&4 to the left of © 
the CG and position 2 at the same distance to the right of the CG. In this 

way Q@ becomes known and together with the values of M and & found previousl 


the necessary information is available to compute @ and Mae 


Since the fluid film bearings are represented by 8 dy..amic coefficients the 
effective bearing stiffness in the horizontal direction differs from the effec! 
stiffness in the vertical direction. Hence, there will be two relative resona) 


associated with each critical speed of the rotor-bearing system. At each of t 
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two resonances the whirl amplitude of the rotor has s relative maximum. The lowest 
resonance is denoted as 'Minor” and the highest resonance is denoted es "Major" 

on the design charts. Which of the two resonances is actually the most pronounced \ 
can be found from the design charts, Figs. D-21 to D-40. Howewer, the speed range | 
between the minor and the major resonances should in general be expected to yield 
relatively large whirl amplitudes such that in the actual application the system 
criticsl speeds will appear as broad amplitude peaks. 


1 
i 
{ 
; 


The design charts are arranged in two groups: Figs. D-l to D-20 giving the critical 
speeds and Figs. D-21 to D-40 giving the corresponding amplitudes. Each chart 
covers a wide variation in the rotor flexibility parameter p and for each value of 

p there are two curves labeled "Minor" and "Major" as explained above. The 


abscissa for the curves is the Sommerfeld number S defined previously. 


In Figs. D-l to D-20 the ordinate gives the value of the speed ratio _ at the 
minor and major resonances associated with the critical speeds. The charts apply 
to both the first and the second critical speeds but such that the value of @ and, 
hence, the values of p and ®, are different depending which critical speed is in- 
vestigated. When p and w are computed the critical speed resonances can be de- 
termined as the intersections between the appropriate design curves and the 
function: 

Oo ge eed 4 

Po DL &? : 


The two corresponding values of the Sommerfeld number can de used directly to 
enter the proper design charts of Figs. D-21 to D-40 and read off the normalized 
rotor amplitude .. x is the rotor amplitude as shown in the Fig. on Page 231. 
Since in an actual application the unbalance is not known (i.e. the magnitude 
of 8 is unknown) the ratio : should be considered as an amplication factor which 
indicates the rotors sensitivity to unbalance. When?p and o. correspond to the 
first critical speed the applied unbalance is a static unbalance whereas for the 


second critical speed the unbalance is a dynamic unbalance. 


When the rotor is very stiff (i.e.p = 0), which is usually the case in gas bear- 
ing machinery, the two first critical speeds (tae rigid body modes, see Page 14) 
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can be calculated from the equation: 


K) x, j 
- | s, : 
| 
1 w{Jl, 2. 22, 4 radians 
On "V2 oot we) ¥ 4 a) mee sec 
where: 
2 yt Ko radians? 
wo = 
t M sec 
t 
2 2 
£2. Kyo, + KY. radians” 
c (I- I,) sec< 
2 
4 (-K,4, > Kt») eadiane” 
mee * MI - IL) sec 
: t P 
and: 
K, »K, Bearing stiffness, lbs/in 
4st, Distances from bearings to CG of rotor, inch 
M. Total mass of rotor, lbs.sec?/in 
I Transverse rotor mass moment of inertia around CG, 1be.in.sec” 
I Polar rotor mass moment of inertia, lba.in.sec” 


(see: Timoshenko: "Vibration Problems in Engineering", D. Van Nostrand Co.N.¥ 


When the supporting bearings are hydrostatic gas bearings tne bearing stiffnesses 


K, and Ky can be found directly from Figs. B-110 to B-122. For hydrodynamic gas 


1 
bearings and hydrodynamic bearings in general the 8 dynamic coefficients can be 
determined from Figs. B-1l to 3-109. The corresponding bearing stiffness to be 
used with the above equation can either be taken in the approximate form: 


~ 1 
K = 2 y + Ky) 
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or, if both the minor and the major resonances are desired: 


I for the minor resonance: K @ Ky 


for the mejor resonance: K Ce Se 


A more exact equation can be expressed as: 


l _ dam -k (ac - aC + H(K oc +K oC 
K= 2 (k,. + Ky + " 


(minus corresponds to the m‘nor resc--nce, plus to the major resonance) 


where A is the positive, real s~ ...on of: 


A * + [Ea eee + KyKyx - dec, = aC, P= at oc of ge 


: (en_- Ky) (at- oc ) + 2K, + yas yy 
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Figure D-1 Critical Speeds of Rotor-Bearing System 
Plain Cylindrical Bearings, L/D = .5, Laminar Film 
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Figure D-2 Critical Speeds of Rotor-Bearing Systen 
Plain Cylindrical Bearings, L/D = 1, Laminar Film 
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Figure D-3 Critical Speeds of Rotor- Bearing System 
4-Axial Groove Bearings, L/D = .5, Laminar Film 
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4-axial Groove Bearings, L/D = 1, Laminar Film 
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Figure D-5 Critical Speeds of Rotor-Bearing System 
Elliptical Bearings, L/D = .5, m= .25, Laminar Film 
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Figure D-6 Critical Speeds of Roter-Bearing System 
Elliptical Bearings, L/D = .5, m= .5, Laminar Film 
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Figure D-7 Critical Speeds of Rotor-Bearing Sy 
Elliptical Bearings, L/D = 1, m= .25, Laminar Film 


245 


~ BOnasearar 
essonenen srainet emomscasmeoratatntt 


WSs 0:3 or.aed* cuss ance 
saenn 7, ‘ care seSssituimniiiisssnae 


mawece-s. 


BSisdiccen ssseess 
ate me pape ehoteme ne tet tf areiar 
sstinimesssetiiins trottitiht a tiinhim 
eaecaseiti uacscscerio auaccersis 
aser amaaessitis 


— 0 a ] 
= ceraee a Ty ~, 9 950001 MG OS Ge 81h181 

es aii: SS esstnieiigne— 2 rere 
oe anan BSS e BaNSal ets: 

@enacess ttt eeeerieei@e 


it fF 


TAT 
ie 
tH [| 
lat at 
Tite 
PHB 
peegeses 
seszerse 
Prete +tts4 
EEz=S 838 
StH 
a 
HA 
3 HE 
saaseg 
es2ss 
Pit +44 
zs2esss 
se eeestie 
wa88 
HEH 
aa 
Hite 
sete 


0.8 meneang Fceetee z Th ite anes seesrasiuiies seasniein 
ees estie ese See 
Siabeabes pases feseiiees ou gate iniy epeeamee sent wives csniin emmcecentt it 

0 7 =osanasas oeoe ses Pett thy fn | ee Ht te 

ss Sta Soar tani eeeeeremeteen ee i 
saass eee ee Sa sscesininegSuisiiiiimsmasasapiiasestaini 
a ss $101} amenene susseneen meen sessss 


As MMOS Bese see me 
it a 
Hot 


L+-4 
3 
atts: Seatiatases 
—- tt +- bende teat 
pane peers ape anes +4 HH 
+ ot 
= PEDES Sa beer tibacs nMNeee @avae ee 
== $2426 thease oe en meae Hitt + tte 
tthe bverthere Heed tt 44 uti m@eaneea 
+t et Pi Hpef pt | meeene 
PPPS S POOR D000 ; Hlt | ag htt eaaees 
ea aeees e8bnestisi Geen ae eceneen 


5 a 


a 
ate 
He 
ane 
wasesonn 
Seles ‘ 5a Oe ne edi: 
eb edoned + + ; ' ‘i 8 
desea ; ha t this tee ge HT 
poe BOCs - see iuteaees ameays 
Re m m3 Senuiesees H 
ipod Gs i DTT aoe esaresiinm ones cesses 
oe + - 
ss Db : Pare re) I O06 a Sosccsss ime eases 
aeccdiimm@ ma mensese 
ft : ode PLt eu 3. ous ’ 
4 - al . te * 


Beet cts or EE ie ORaes teens FEE 
sszsiith 
Seta EEE Hence 
= ae 
0. see ies ae EE 
= ; be bane Ee 


BSS (a erthe 
fa : 
a ree ‘Rie TC 


et 
1 


==... INAS Ri sessel 
(J ages 
amma asscessisuie 
t jot ttt tty 
manesas ae 
jeeurs 


=aee HHH etieS 


2 4 6 *1.0 


Figure D-8 Critical Speeds of Rotor-Bearing System 
Elliptical Bearings, L/D = 1, m6 .5, Laminar Film 
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Figure D-9 Critical Speeds of Rotor-Bearing System 
6 Shoe Tilting Pad Bearings, L/D m= .5, Load between Pads 
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Figure D-10 Critical Speeds of Rotor-Bearing Systes 
5 Shoe Tilting Pad Bearings, L/D = .5, Load between Pads 
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Figure D-11 Critical Speeds of Rotor-Bearing System 
& Shoe Tilting Pad Bearings, L/D = .75, Load between Pads 
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_ 6 Shoe Tilting Pad Bearings, L/D = .5, Load on Pad { 
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Figure D-13 Critical Speeds of Rotor-Bearing System 
5 Shoe Tilting Pad Bearings, L/D = .5, Load on Pad 
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Figure D-14 Critical Speeds of kotor-Bearing System 
4 Shoe Tilting 2ad Bearings, L/D = .75, Load on Pad | 
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Figure D-i5 Critical Speeds of Rotor-Bearing Systes 
100 Degree Partial Bearings, L/D = 1, Turbulent Pils 
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Figure D-16 Critical Speeds of Rotor-Bearing System 
160 Degree Partial Bearings, L/Del, Turbulent Fila 
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Figure D-17 Critical Speeds of Rotor-Bearing Syetem 
100 Degree Partial Bearings, L/D=l, Turbulent Film 
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Figure D-18 Critical Speeds of Rotor-Bearing System 
Plain Cylindrical Bearings. L/D = 1, Turbulent Film 
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Figure D-19 Critical Speeds of Rotor-Bearing System 
Plain Cylindrical Bearings, L/D = 1, Turbulent Film 
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D-20 Critical Speeds of Rotor-Bearing System 


Figure 
Plain Cylindcical Bearings, L/D = 1, Turbulent Film 
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Figure D-21 Dimensionless Rotor Amplitude at the Critical Speed 
Plain Cylindrical Bearings, L/D=.5, Laminar Film 
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Figure D-22 Dimensionless Rotor Amplitude at the Critical Speed 
Plain Cylindrical Bearings, L/D = 1, Laminar Film 
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Figure D-23 Dimensionless Rotor Amplitude at the Critical Speed 
4-Axial Groove Bearings, L/D = .5, Laminar Film 
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Figure D-25 Dimensionless Rotor Amplitude at the Critical Speed 
Elliptical Bearings, L/D = .5, m= .25, Laminar Film 
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Figure D-26 Dimensionless Rotor Amplitude at the Critical Speed 
Ellipticai Bearings, L/D = .5, m= .5, Laminar File 
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Figure D-27 Dimensionless Rotor Amplitude at the Critical Speed 
Elliptical Bearings, L/D= 1, m= .25, Laminar Fila 
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Figurs D-28 Dimensionless Rotor Amplitude at the Critical Speed 
Elliptical Bearings, L/D = 1, m= .5, Laminar Fila 
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Figure D-29 Dimensionless Rotor Amplitude at the Critical Speed 
6 Shoe Tilting Pad Bearings, L/D = .5, Load between Pads 
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Figure D-30 Dimensionless Rotor Amplitude at the Critical Speed 
5 Shoe Tilting Pad Bearings, L/D = .5, Load between Pads 
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Figure D-31 Dimensionless Rotor Amplitude at the Critical Speed 
4 Shoe Tilting Pad Bearings, L/D = .75, a between Pads 
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Figure D-32 Dimensionless Rotor Amplitude at the Critical Speed 
6 Shoe Tilting Pad Bearings, L/D ~ .5, Load cn Pad 
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Figure D-33 Dimensionless Rotor Amplitude at the Critical Speed 
5 Shoe Tilting Pad Bearings, L/D = .5, Load on Pad. 
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Figure D-34 Dimensionless Rotor Amplitude at the Critical Speed 
4 Shoe Tilting Pad Bearings, L/D = .75, Load on Pad 
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Figure D-35 Dimensionless Rotor Amplitude at the Critical Speed 
100 Degree Partial Bearings, L/D = 1, Turbulent Film [ : 
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Figure D-36 Dimensionless Rotur Amplitude at the Critical Speed 
100 Degree Partial Bearings, L/D = 1, Turbulent Film 
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Figure D-37 Dimensionless Rotor Amplitude at the Critical Speed 
100 Degree Partial Bearings, L/D = 1, Turbulent Film 
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Figure D-38 Dimensionless Rotor Amplitude at the Critical Speed 
Plain Cylindrical Bearings, L/D = 1, Turbulent Film 
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i Figure D-39 Dimensionless Rotor Amplitude at the Critical Speed 
Plain Cylindrical Bearings, L/D = 1, Turbulent Film 
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Figure D-40 Dimensionless Rotor Amplitude at the Critical Speed 
Plain Cylindrical Bearings, L/D = 1, Turbulent File 
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NUMERICAL EXAMPLES, 


Three numerical examples are given to illustrate the use of the design charts. 
For this purpose three different rotor-bearing systems have been selected: one 
configuration where the rotor is flexible and the bearings are oil lubricated, 
and two configurations with a rigid rotor where the bearings in one case operate 
with a turbulent film and in the second case the bearings are gas lubricated. 


The configurations are chosen to be typical rather than representing any parti- 
cular application. The calculations are concerned with the dynamics of the 
systems: the dynamic bearing coefficients, the critical speeds and the stability. 
The steady-state performance is in general not considered. 

1 otor m r l 


An electrical motor is mounted on a shaft between two bearings. The motor drives 


a pump wheel which is overhung at one end: 


The following data are given: — 1 
Total rotor weight: 155 lbs. 
Bearing reactions: 60 lbs.and 95 Ibs. 
Bearing span: 28 inch 
Shaft diameter: 1 inch 
Youngs modulus for shaft: 3 - 107 Lbs/in” 
5,000 RPM. 


Design speed: 


Journal bearing type: 4-axial groove bearings 


Bearing diameter: D=1 inch 
Bearing length: L = 0.5 inch 
280 
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Radial clearance: Cc = 0.001 inch 

Oil: SAE 10 at 120°F 

O11 viscosity: = 24 centipoise = 3.5 - 10° Lbs. sec/in? 
Oil density: ep = 0.0307 Lbs/in? = 7.95 °° 1079 lbe.sec?/in’ 


It is desired to investigate the critical speeds, the unbalance response and the 
stability of the given rotor-bearing system to determine if the performance of the 


system is satisfactory with the selected rotor and bearing dimensions. 


For analysis purposes the rotor car. be considered to consist of two masses: the 
pump wheel with a weight of 30 lbs. and the motor rotor, weighing 125 lbs. with 
a transverse mass moment of inertia of 6,940 lbs. in® and a polar moment of inertia 


of 260 lbs.in?. The shaft is 1 inch in diameter with dimensions as shown in the 
figure. 


Consider first the stiffness of the rotor. The rotor has three degrees-of-freedos 
namely the amplitude x) of the motor rotor, the amplitude Xo of the pump wheel 
and the angular displacement 9 of the motor rotor. Let the three quantities be 
identified by subscript 1, 2 and 3, respectively, and let @ denote an influence 
coefficient such that a, means the shaft deflection at the wheel due to a 1 1b. 
force at the same position. Similarly, a5, means the deflection at the wheel due 
to a Ll lb.in. moment at the CG of the motor rotor. From standard beam formulas 
(See: R.J. Roark: "Formulas for Stress and Strain", McGraw-dill Book Co., New Yor 


the values of the influence coefficients are determined to be: 


“6 in -8.15- 107° 2 ” <1 
@,57.24- 10° FS, a, = -8.15- 10° >, a= 0 lbs 


= 0.0223 - 10°° 7adians 


© 1976 de aoe a0? be 
a = 28.05 10 ? a 0.771 10 lbs, 43 lbs.in. 


22 lbs 


7 


The values are basedonaYoung's modulus of 3 - 10 lbs/in” and a cross-sectional 


shaft moment of inertia of 0.0491 in’, corresponding to a l inch diameter. The 
three inertias become: . 
mass of motor rotor = 125/386 = 0.324 Lbe.sec /in 


mass of pump wheel = 30/386 = 0.0777 lbs.sec2/in 2 
mass moment of inertia of motor rotor = (6,940-270)/386 = 17.3 lbs.in.sec 
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The equations of motion become: 


x, 7 1: 0.324 . x * A>, 0.0777 . Ro -A3- 17.3 @ 
j 


Ry @ Ao: 0.324. %, 7 Boo 0.0777 . Ry 7 Ay + 17.3 6 


@ = -a,,.- 0.324 . x, - a), 0.0777 . x) - a, . 17.36 
Let the rotor speed be w Eadiens 


where * means the x- acceleration, ets. 
whereby x" ox, etc., and subetitute the previously derived values for the 


influence coefficients to get: 
= (0 


(2.345 . aw - 10°) x, - 0.633. a. x, 


“2.64 . wo x, + (2.18 . aw - 10°) x - 13.33.07 @=0 


-0.0599 . wx, + (0.386 . a - 10°) 9 =0 


From the last equation: 
0.0599 w 
©" 0.386 . a - 10 
in which case the determinant for the Xy7 Ro equations can be written: 
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10° 
(2.345 - 12 -0.633 
~2.64 (2.18 - 10° ~ 23:32..0.0599)) 2 9 
2 (0.386 = —) 
a ries w* 


By trial and error the two lowest solutions for @ become: 


w = 524 44 5 5,000 BPM 
aec 


lst rotor resonance: 


rad 
w= 865 pace 8,250 RPM 


2nd rotor resonance: 
These resonances will be lowered due to the effect of the flexibility of the 
bearing film. To get an estimate of this effect it may be noted that the first 
of the two rotor resonances basically corresponds to the resonance of the motor 
rotor on the shaft. Hence, with a mass of 0,324 lbe. sec” /in. the total rotor 
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stiffness yielding a frequency of 524 isdians/sec becomes 86,000 lbs/in. Since 
the design charts consider the rotor to be symnetric with half the rotor mass 
assigned to each of the two bearings the rotor stiffness to be used in the charts 


should be taken as half of the total stiffness value, 1.e.: 


lba 
K.. 2 43,000 in 


The effective influence coefficient a (see page231) becomes: 


5 
‘ 


1 -5 in 
a K = 2.325 . 10 ibe 


and the rotor mass per bearing is 0.162 lbs. sec? /in. 


The bearings are oil lubricated and are of the 4-axial groove type. The journal 
diameter is: D = 1 inch and the bearing length is: L = 0.5 in. The radial 
Clearance is: C = 0.001 inch and the oil is SAE 10 at 120 F with: 
Viscosity: wp = 24 centipoise = 3.5. 107° Lbe.sec/in® 
3 <5 2,,..4 
Density: p = 0.0307 lbs/in’ « 7.95 . 10  lba.sec /in 


With a bearing load of: W = 62.5 lbs the Sommerfeld number becomes: 


-6 2 
3.5.10 .~N.1. 0.5 ,0.5- -3 
62.5 G00 ee 
To determine the first critical speed use Fig. D-3. With a= 524 kad 


the relationship between the Sommerfeld Number and the speed ratio is: 


8 = 0.585 oe) 
n 


The rotor flexibility parameter is: 
-5 


_ 62.5, 2.325. 10> _ 
e 0.001 1.45 


Plot the curve of $ versus a, Fig. D-3 and determine the intersections with 


the curves corresponding to the calculated P-value: 


for major resonance: 5 = 0.98 , § = 0.573 
n 


for minor resonance: = = 0.88 , § = 0.515 
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The corresponding response is found from Pig. D-23 as: 
for major resonance: 4 =7 

x 

& 


for minor resonance: 


Hence the minor resonance is the most pronounced and the first critical speed 


of the system becomes (in approximation): 
let critical speed = 5,000 . 0.88 = 4,400 RPM 


is predominantly controlled by the overhung pump wheel 


The second critical speed 
However, it is close 


and actually the design charts do not apply to thie case. 
enough to the first critical speed that the same reduction factor can be used 


so that: 
2nd critical speed = 8,250 . 0.88 = 7,300 RPM 


The above calculations are only approximate and should be checked by more elabor- 


ate methods such as a computer program. For this purpose it is necessary to 


know the dynamic bearing coefficients. From the figure on page280 the bearing 
reactions are found to be: 


for left hand bearing: W= 60 lbs. i.e., 2 = 6. 10° 1bs/in 


= 9.5 . 10° lbs/in 


Oe als 


for right hand bearing: W = 95 lbs. i.e., 
At 4,400 RPM the Sommerfeld Number becomes: 
for left hand bearing: S = 0.535 


for right hand bearing: S = 0.338 


Using Figs. B-5 and B-6 the dynamic coefficients can be computed: 


for left hand bearing 


CU ae Slbs y . Ibs 
K m 2.46 . 10 in? Ky 1.8. 10 in” Kx , 2,400 in” Ky 


= ibe.sec a ibs.sec 
4.2. tn Cy 10.4 in 


1 


= 5.8. 10° ob 


ibs.sec = -195 Lbs.sec 
C= 430 SHEE, ct 7195 ST, 
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for right hand bearing 


5 lbg 5 lbs lbs 
K. 7 4-65 . 10 in” re 3.09 . 10 K = 16,000 7, 


in’ yx 
5 lbs 
Ky = 1.06 . 10 i. 


lbs.sec Pet lbs.sec “3 lbs.sec " lbs. sec 
oS e 595 ta P oy 340 ante a w «5.2 a C, 2.1 in 


These values can be used directly as input to a computer program for calculating 
the detailed unbalance response of the rotor-bearing system and to check the 
stability at the selected speed. If the response or the stability is investi- 


gated at other speeds the above procedure can be repeated to obtain the value of 


the bearing coefficients at the new speed. 


The stability of the rotor-bearing system can also be investigated by meens of 

the present design charts. Since it is found above that the lowest resonance 

(the first critical speed) is basically governed by the motor rotor mass, the 

mass per bearing is set equal to: M = 0.162 lbs.sec’ /in and the rotor stiffness | 
is K,= 43,000 ibe The corresponding bearing load is; W = 62.5 lbs. Compute 


the dimensionless rotor mass parameter: 


V CMW Vo.001 . 0.162 . 62.5 = 0.23 


ubL @)? 3.5.10°° . 1.0.5. (0.5. 10)" 


From Fig. C-2 the corresponding Sommerfeld Number is: S$ = 0.94 which yields a 
speed of: N = 134 RPS = 8,040 RPM. This is the speed ct the threshold of 
instability if the rotor is rigid. However, since the rotor is flexible the 
actual threshold speed will be lower. Following the procedure outlined on page a 
203 the reduction due to the rotor flexibility may be computed as given below: , 


cK 
- Q 1 
k= Ew 2.001 43,000 6 9 69g 


from Pig. C-2: 


N,; RPS re} M K = SA, 7 M, -Vik/ (ke + K,) 
100 0.7 0.355 0.6085 0.259 
117 0.82 0.28 0.521 0.211 
134 0.94 0.23 0.461 0.178 
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The last column gives the modified dimensionless rotor mase parameter as 4 
: function of Sommerfeld number. With the above calculated mass parameter of 0.23 | 
: the Sommerfeld number at the threshold of instability becomes: § = 0.78 which 


; 
corresponds to a speed of: 


Returning to the previously calculated response at the first critical speed: 
xe 23 an estimate can be made of the level of balance required for the motor 


& 


| 
Speed at threshold of instability = 111 RPS = 6,650 RPM | 
rotor. Let the maximum tolerated rotor displacement be: x = 0.0005 in. The | 


eccentricity of the motor rotor is 5 and may be expressed as: 


g » Wnbalance, lbs.in) | Unbalance ‘ . . | 


(Rotor weight, lbs) 125 


Thus: | | 


Max. unbalance = 125 .—%& = 125. 2S = 0.0027 Lbs.in = 0.043 oz. inch | 


: 
The required balance level of the pump wheel cannot be analyzed by means of the 
design charts. For this purpose a computer program may be applied. However, the | 
calculations performed above indicate that the maxima unbalance of the pump 
wheel should not exceed approximately 0.01 oz.in. | 


Finally, the operating speed of the unit is 5,000 RPM. Computing the corresponding | 
Sommerfeld number the bearing eccentricity ratio can be determined from Fig. A-1 
and the triction power loss can be estimated from the equation given on page 24 : 


| 
Ss € Min. Film Thicknesa, inch Power Loss, HP | 
Left hand bearing 0.668 0.555 0.900445 0.069 
Right hand bearing 0.384 0.625 0.000375 0,073 
Total 0.14 HP 


The major conclusion to de derived from the review of the dynamic performance 

of the unit is that the shaft is too flexible. The operating speed is too close : 
to the first critical speed and the bearings are sufficiently stiff that any. | | 
modification of them will not contribute to changing the critical speed. The 

very flexible shaft also makes the rotor quite sensitive to any unbalance. _ | 
Hence, the design of the rotor mst be changed. The most effective way is to 
increase the shaft diameter if the resulting increase in the power loss is 
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b. Rigid Rotor, Turbulent Bearing Fila 


A emall rigid rotor is supported in two Mercury lubricated bearings. The 
system has the following data: 


Total rotor weight: 10 lbs. 

Transverse minus polar mass woment of inertia of rotor: (I-1,) = age 3 lbs. in®= 
0.23 lbs.sec in. 

Bearing reactions: 5 iba. and 5 lbs. 

Bearing span: 10 inch 

Design speed: 40,000 RPM 

Journal bearing diameter: D = 0.625 in. 

Bearing length: L = 0.625 in. 

Radial clearance: C = 0.0007 in. 

Lubricant: Mercury 

Lubricant temperature: 400 F 

Lubricant viscosity: ,: = 1.51. 107” Lbs. sec/in® 

Lubricant density: Pp = 0.475 1be/in? = 1.23 . 107? 1be.sec*/in* 


It is derired to select the bearing type for the application. 


The rotor weighs 10 lbs and can be taken as being symmetric such that the load 
on each bearing is 5 lbs and the rotor mass per bearing becomes: 


M = = .013 lbs. sec? /in 


ae 
386 
The journal diameter is: D = 0.625 in. and the bearing length is: L = 0.625. The 
radial clearance is: C = 0.0007 in. which is relatively large but needed to 
accommodate thermal expansions and to provide adequate cooling. With these 
dimensions and the known lubricant properties the bearing parameters can be 
calculated. Denoting the speed as N (in RPS) we get: 


-7 . 
1.5 . 10 0.625 . 0.625 3125.2 <3 
. eS ee 
Sommerfeld pumber: & 5 ese m2.46 . 10 - N 


-3 
sajuates miapees: - see wi 222 s HOT Mag 02625400007 a1. a 


The operating speed is 40,000 RPM. Investigate first if it is possible to use 
a simple bearing type like a plain cylindrical bearing or a grooved bearing. 
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For study purposes the latter type can be represented by a 100 degree partial 


bearing. 


It seems evident that the most difficult problem is instability. Hence, calcu~ 


late the dimensionless rotor mass parameter: 


Vow Vo.0007 . 0.013 . 5 wane 
“R27 1-51. 10-7. 0.625. 0.625, 0.3105,2 "9 
uD.) 0007? 


{ 
. 


Entering Fig. C-9 with this value and using the curve for a 100 degree partial 


bearing gives: S$ = 0.39. The corresponding speed becomes: 


0.39 
Speed at instability threshold = 0.00246 " 159 RPS = 9.500 RPM 


The corresponding Reynolds number is 1770 which is sufficiently close to the 


value of 1663 for which Fig. C-9 is valid. Thus, a 100 degree partial bearing 


yields a threshold speed of 9,500 RPM which is far too low for the application. 


Next, enter Fig. C-1l1 with the calculated rotor mass parameters an? use the 


curve for the plain cylindrical bearing to get: 


0.64 ji 
$= 0.64, i.e. N= “7to7e = 260 RPS = 15,600 REM 


The corresponding Reynolds number becomes 2910 which agrees well with the value 


of 3326 for which Fig. C-1l is valid. Hence, the threshold speed is 15,600 RPM 


for a plain cylindrical bearing which is not acceptable. 


It can, therefore, be concluded that no "fixed geometry" bearing is likely to 
ensure stable operation. Among the other possible bearing types (the floating 


ring bearing, the externally pressurized bearing, etc.) the tilting pad bearing 


is the one for which most experience is available. It wiil, therefore, be 


investigated if the tilting pad bearing is adequate for the present application. 
To check its stability it is sufficient to check if the pivoted shoes follow the 
For this purpose Fig. 3-45 can be employed to get an estimate 
For a speed of 40,000: RPM = 
The corresponding critical 


journal motion. 
even if the chart assumes the film to be laminar. 
667 RPS the Sommerfeld number becomes: S$ = 1.64. 


mass is found to be: 
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or 
Mie 7 0-367 Ibs.sec*/in = mace inertia of shoe 
: (0.3125) 
Hence, 


Critical masa moment of inertia of shoe = 0.0359 liwsitaceec’ = 13.85 lbs 


If the shoe is a 90 degree arc its "curved" length is: 7]. 0.625/4 = 0.49 in. 
With a length of 0.625 inch, a thickness of 0.25 inch and a density of 


0.283 lbs/in? its mass moment of inertia is approximately: 


2 
Mass moment of inertia of ehoe = 0.49 . 0.625 . 0.25 . 0.283 . o.%) = 


0.000435 Ibs. in’ 


a 
which is far below the critical value. Hence, the tilting pad begrings will | 


stable. 


The dynamic bearing coefficients can be found from Figs. 8-80 to 8-86. Once 
coefficients are determined the critical speeds may be cal.ulated directly sii 
the tilting pad bearing has no cross-coupling terms. The procedure is illus« 


trated below. It consisty of the following steps: 


1. assume a rotor speed N, RPS 

2. calculate the corresponding Reynolds number and Sommerfeld number 

3. enter the appropriate design chart (or interpolate) to find the 
spring coefficients 

4. calculate the corresponding critical speed 

5. if the calculated critical speed differs greatly from the assumed 


speed, repeat the procedure. 
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The rotor mass per bearing {80.013 lbe.sec7/ia, the transverse wees woment of 


inertia minus the polar mass moment of inertia of the rotor is: (I-1,) = 88.8 
ibe ane =0.23 lbs. in. sec” and the bearing span is: #4 = 10 in. If the bearing 


atiffness is called K ibs the critical speeds become: i 


1 
lst critical speed = 2x 013 
1 


= 2 
: iL | 2ns" L JK 
2nd critica. speed = x 0.23 Sag 0046 BPS 


The bearing reaction is: W= 5 lbs. and the radial clearance is: C = .0007 in. 


such that: 


3 lbs 


ys ; 
i Oe Le 


Following the outlined procedure the calculations yield: 
N= 125 2PS S$ =0,308 Re = 1400 


ere lbs lbs, gec 
from Fig. B-81: K = 6,400 7, CLs 21.8 ia , 


i.e. lst critical speed, minor resonance = i »f5 400 ws 112 BPS = 6,700 RE 


N = 180 RPS S =0.443 Re = 2020 
lbs : lbs.sec 
from Figs. B-81 and B-82: nyo 15,350 7 , c 55.6 


i.e. Ist critical speed, major resonance = 7 ‘ 15.390 = 173 RPS = 10,400 ! 


lbs lbs .sec 
from Figs. B-81 and B-82: ad 5,200 in? C2 18.6 are | 
i.e. 2nd critical speed, minor resonance = - 2.200 = 169 RPS = 10,200 8 
N.=_250 RPS S =Q615 Re = 2800 
s Aa m lbs ~ lbs .sec 
from Figs. B-82 and B-83: mS 10,930 jan 2 Cy 52.7 in 


i.e. 2nd critical speed, major resonance = i = f'Sa230 = 245 BPS = 14,704 


{ 


It should be noted that the bearing damping is sufficiently high that the systé 
is critically damped st all the calculated resonances. Thus, at the first crit 
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speed: | 18.2 ibe ses (ist criticsel, minor) 
Critical damping = 2 ‘¥K° .013 = HSecaee 
28.2 ae (1st critical, major) 


Similarly for the second critical speed. Thus, the critical speeds will have | 
influence on the rotor performance. Lc can,therefore, ba concluded that the ¢ 
pad bearing ensures satisfactory performance of the system from the point of v 
of rotor-bearing dynamics. 


c. Rigid Rotor. Hydrostatic Gas Bearings 


A rotor has the following data: 


weight: 90 ibs. 

transverse mass moment of inertia at CG: 7,230 ibs.in® 
tearing span = 20.7 in. 

distance from bearing #1 to CG = 11.5 in 

journal diameter = 2.5 in 


max. journal bearing length = 3.75 in. 


The operating conditions are: 


normal operating speed = 32.000 RPM = 533 cps 
max. overspeed = 36,000 RPM = 600 cps 

max. available supply pressure = 325 psia 
gas = air at 120°F 


It is desired to select the dimensions of the hydrostatic bearings for this 


application. The bearing reactions are found to be 


—= eo 


Brg. #1 11.5 9.2 Brg. #2 


90 lbs. 
for bearing 1: W = 40 lbs. 


for bearing 2: W = 50 Ibs. 


Hence, use W = 50 lbs. for design purposes. Since the rotor is heavy for the 
required top speed it is necessary tc have the bearing as stiff as possible & 
order to avoid hydrodynamic instability. Therefore, the journal bearings are 
chosen to have double plane admission. From the approximate relationships: 
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With L = 3.75 in. it 1s evident that the required load carrying capacity is easily 


; achieved. Turning to the bearing stiffness it is necessary to have: 


lowest critical speed > 0.6 - (max. speed) 


The lowest critical speed can be estimated by assigning a rotor mass equal to ; 


the load of the heaviest loaded bearing, i.e. bearing No. 2. Hence, the re- 


quired stiffness K may be estimated from, 


: L (— : 
ox ¥50/386.07. 7 70 S00 = 280 eps | 
K> 6.6 10° ibe 


Choose the optimum value of the restrictor ratio: ¢ 
Ag =0.7 


and let the bearing have the maximum allowable length: 


; L = 3.75 in. 

: i.e. L. 3.75 
D257 ! 
t L 
and b= of au7s 


stif fneq 


For an estimated operating pressure ratio of: p/P — 2 the dimensionless 
is determined by interpolation between Figs. B-119 and B-122: 


i 1+ oo? PPD | 
i 


The journal bearing clearance is selected such that: 
\ 

gc. - 1072 
R 0.6 - 10 


ia C 90.6 + 1072 + 1.25 =0.75 + 10") in, 
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’ becomes: 


To get the maximum stiffness the feeder holes should be orifice restricted with 


- ae little inherent compensation as practical, i.e. 5« 0. Then the actual 


stiffness becomes: 


1 
e 3 - 
R975 107) (P ,-P,) 3.75 + 2.5 +0.47 = 5.88 + 19 , Pr, 


Therefore 


5 
6.6 * 10 
Pa Po pd mad pst 


or: 
P, > 113 + 120 = 233 psia 


To allow some margin set the actual operating supply pressure equal to 275 psie. 
Since the maxivum available supply pressure is 325 psia the required stiffness 
can be obtained with the chosen bearing dimensions. 


The lubricant is air at 120°F with the properties: 
-9 1bs.sec 


in” ‘ 
in2/sec2 °R 


viscosity: p= 2.8 ° 10 


gas constant: R= 2.472 > 10° 
2 
ie. QT = 2.472 + 107+ (460 + 120) = 1.434 + 10° AR 
sec 
or VQ@T@= 1.198 - 10° Sh 

sec 
To select the number of feeding holes N and the orifice radius a for the joumal 

bearing use the definition of the restrictor coefficient A: 
2 Pc 


a a A 4 


148 ouvert = 


For N= 16 and & = 0, a= 5.8 +10 in. Since it is not feasible to have 8 


exactly equal to zero, set: 


‘O75 ° §.3823-10. 


a= ,006 in (i.e. orifice diameter #0.012 in) 


with N= 16 


With a feeder hole diameter d =0.125 in. the inherent compensation factor 


2 2 
me 0.006 < 
B= Gc 70.125 -0.75 1073 70-384 
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x 
so that the reduction in stiffness is lese than 5 percent. Check the value of } 
the reat?ctor coefficient: 


-9 2 and 
Ag = S28 + 10°: 16 + (0.006)" 1.198 - 10" 9 96 5 79 
t 275 -(0.75 - 1073) V1 + .3847 


which agrees with the selected value. 


Let the distance from the orifice to the bearing surface be0.050 in. Then with 
a hole diameterof 0.125 in. the feeder hole volume becomes: ve 0.613 - 10°? in? 


and the feeder hole volume ratio is 


We ft 16 + 4613+ 107? mare 
aDLC / -3 : 
x /2.5 + 3.75 + .75 + 10 


which is toohigh by an order of magnitude from the point of view of susceptibility 
to pneumatic hammer. Instead make the feeder hole diameter: 
#9,035 i 
and provide the feeder hole with a recess with a diameter of .125 inch and a 
depth of0.002 inch. 


The feeder volume ratio then becomes: 


( 
xDLC =0.0525 


which is acceptable. Yroviding the recess also changes the inherent compensation 


factor which becomes: 


& =0.536 


Based on the already obtained value for the dimensionless stiffness, the actual 
stiffness becomes: 


2 i 

144 + (0.536) : 

ny iba 
= =3 ° (275- - 3. ° 2.5 °0.47 = 8.43°10 


\ 
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Assume a support stiffness of 3° 10° ree. so that the combined stiffness of | 


the gas film and the support becomes: 


combined stiffness = 6.58 - 10° ibe 


To check the critical speeds of the rotor assume the rotor to be rigid. The masd 


and mass moment of inertia of the rotor are computed from the given data: ; 
i 


2 
= lba ,8ec i 
0,233 in 


i 
i 


; pee Lt 
rotor mass: x." 386.07 


(21230-8660), 16 5 ibs.in. sec” 


rotor mass moment of inertia: (I-1,) = 386.07 


Using the equations on Page 235 with kK °k,° 6.58 °* 10° pba and t, «= 11.5 in 


and s, = 9.2 in., we get: 


. 5 2 
we a 26:58 10" 5 gg. 196 zad 


t 0.233 sec? 
5 2 2 
we a $298. 2 10" (115° + 9.2%) _ 8.65 - 10° ad 
c 16.5 2 
sec 
12 rad¢ 


a = (6:58 + 10°(-11.5 + 9.2) )? =0.596 + 10 


ct 16.5 °6.233 sec 


and the two lowest critical speeds become: 


w= 10° - , (8.65 + 5.65) +y} (8.65 - 5.65)7 +0.596 


372 cps (translatory mode) 


Rigid body critical speeds= , 
472 cps (conical mode) 


Hence: 


The lowest critical speed (i.e. 372 cps) is higher than the required 60 per cent 
of the maximum rotor speed. The second critical speed is 13 per cent below the 
normal operating speed of 533 cps and should, therefore, not seriously interfere 


with the rotor operation. 


To get an idea of the response at the critical speed determine the damping co- 
efficient. The dimensionless value is obtained by interpolation between Figs. 
B-120, B-121, B-123 and B-124. 
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L Wy Dimensionless 


1 0 4.0 
lL 0.1 3.7 
2 0 12.0 
Zs Ou 11.8 | 
1.5 0.05 7.9 , | 


where the last line is determined by linear interpolation. The actual dasaping 


B= 28° 10% - 3.75 ¢-t28—)? - 7.90383 ersee 


becomes: 


The lowest critical speed (the translatory mode) is closely determined by 


7 my bse ee Jo2  tadtane 
n M t sec 
2 


(i.e. we = 5.65 > 10° a or a = 385 cps which is very close to the calculated 


372 cps for the translatory mode). 


3 2 
sonance of a mass0. 233 Sbe.s8e- supported by a spring of: 2 - 6.58 * 10 


10° ibs The critical damping is defined as: 


F - 10° .« lbs.sec 
Biritical™ 2 1.316 10° -0.233 = 1106 “T+ 


Thus, the log-decrement 5 is approximately: 


Bw SB a7 


Boritical 


Hence, this mode can be considered as the re- 
3 2 1.316 « 


and the magnification factor (the Q-factor) is closely given by: 


Q-factor = = = 1.45 


Note: the Q-factor gives the ratio between the amplitude at resonance and the 


amplitude if there was no resonance, i.e. it measures the amplification due to 
the resonance such that a Q-factor of 1 means no amplification. ‘hus, the rotor 


is rather heavily damped. 
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The operating eccentricity ratio can be calculated from: 


és = =0.079 


Finally, a summary of the bearing dimensions and the most important of the design 


data shall be given: 


Resting Data 


Diameter: D = 2,5 in. 

Length: L = 3.75 in. 

Length between admission planes: L, » 1.875 in. 

Radial clearance: C #0.00075 in. 

Number of admission planes: 2 

Number of feeder holes: N= 16 (8 holes per admisecion plane) 

Orifice diamter: 2a =0.012 in. 

Feeder hole diameter: d =0.035 in. 

Recess diameter: =Q.125 in 

Recess depth: #=0.002 in. 

Length of feeder hole: * 20.050 in. 

Supply pressure: ?. = 275 psia 

Radial stiffness: K = 840,000 lbs/in. 

Operating eccentricity ratio: € =0.079 

Critical speeds of rotor: 372 cps = 22,300 RPM (translatory mode) 
472 cps = 28,300 RPM (conical mode) 
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